1/2 


AD-A173  342 
UNCLASSIFIED 


PROCESSING  SCIENCE  TO  INCREASE  THE  RELIABILITY  OF 
CERAHICS(U)  ROCKWELL  INTERNATIONAL  THOUSAND  OAKS  CA 
SCIENCE  CENTER  B  J  KELLETT  ET  AL.  SEP  8C  SC41S.  FR 
N00014-84-C-8298  F/G  11/2  NL 


1,0  S  la  itt 


II  S."*  HS 


m  n  1.6 


MICROCOPY  RESOLUTION  TEST  CHART 
NATIONAL  BUREAU  OF  STANDARDS- 1963-A 


NIC  ^ 


AD- A 173  342 


SCS410.I 
Copy  No. 


PROCESSING  SCIENCE  TO  INCREASE 
THE  RELIABILITY  OF  CERAMICS 

TECHNICAL  REPORT  NO.  1 

THERMODYNAMICS  OF  DENSIFICATION,  PART  I:  SINTERING  OF  SIMPLE  PARTICLE 
ARRAYS,  EQUILIBRIUM  CONFIGURATIONS,  PORE  STABILITY,  AND  SHRINKAGE 

B.J.  Kellett  and  F.F.  Lange 
TECHNICAL  REPORT  NO.  2 

THERMODYNAMICS  OF  DENSIFICATION,  PART  II:  GRAIN  GROWTH 
IN  POROUS  COMPACTS  AND  RELATION  TO  DENSIFICATION 

F.F.  Lange  and  B.J.  Kellett 
TECHNICAL  REPORT  NO.  3 

THERMODYNAMICS  OF  DENSIFICATION,  PART  III:  EXPERIMENTAL 
RELATION  BETWEEN  GRAIN  GROWTH  AND  PORE  CLOSURE 

B.J.  Kellett  and  F.F.  Lange 
TECHNICAL  REPORT  NO.  4 

EXPERIMENTS  ON  PORE  CLOSURE  DURING  HIPing  AND  FORGING 
B.J.  Kellett  and  F.F.  Lange 


FINAL  REPORT  FOR  THE  PERIOD 
July  1,  1984  through  March  31,  1986 


CONTRACT  NO.  N00014-84-C-0298 
PROJECT  NO.  4310037 


Prepared  for 

Office  of  Naval  Reaearch 
800  North  Quincy  Street 
Arlington,  Virginia  22217 

F.  Lange  and  B.  Kellett 
Principal  Investigators 

SEPTEMBER  1986 


s 

t- 


DTIC 

EL.ECTE 
OCT  2  2  686 


Approved  for  public  release;  distribution  unlimited 


4* 


Rockwell  International 


Science  Center 


FTOP 


/ 


TABLE  OF  CONTENTS  ' 


Rockwell  International 

Science  Center 


SC5410.FR 


Page 


Technical  Report  No.  1 

THERMODYNAMICS  OF  DENS I F I  CAT I  ON ,  PART  I:  SINTERING  OF 
SIMPLE  PARTICLE  ARRAYS,  EQUILIBRIUM  CONFIGURATIONS, 

PORE  STABILITY,  AND  SHRINKAGE. .  1 


Technical  Report  No.  2 

THERMODYNAMICS  OF  DENS I F I  CAT I ON ,  PART  II:  GRAIN  GROWTH  IN 
POROUS  COMPACTS  AND  RELATION  TO  DENSIFICATION : .  43 


Technical  Report  No.  3 

THERMODYNAMICS  OF  DENSIFICATION,  PART  III:  EXPERIMENTAL 
RELATION  BETWEEN  GRAIN  GROWTH  AND  PORE  CLOSURE .  79 


iXti 


v  -A 


Technical  Report  No.  4 

EXPERIMENTS  ON  PORE  CLOSURE  DURING  HIPing  AND  FORGING 


101 


SC5410.FR 


THERMODYNAMICS  OF  DENSIFICATION,  PART  I:  SINTERING  OF  SIMPLE 
PARTICLE  ARRAYS,  EQUILIBRIUM  CONFIGURATIONS, 

PORE  STABILITY,  AND  SHRINKAGE 

B.  Kellett  and  F.F.  Lange 

Department  of  Materials  Science  and  Engineering 
University  of  California,  Los  Angeles 

and 

Structural  Ceramics  Department 
Rockwell  International  Science  Center 


Rockwell  International 

Science  Center 


THERMODYNAMICS  OF  DENSIFICATION,  PART  I:  SINTERING  OF  SIMPLE 
PARTICLE  ARRAYS,  EQUILIBRIUM  CONFIGURATIONS, 

PORE  STABILITY,  AND  SHRINKAGE 

B.  Kellett  and  F.  F.  Lange 


Department  of  Materials  Science  and  Engineering 
University  of  California,  Los  Angeles 


Structural  Ceramics  Department 
Rockwell  International  Science  Center 


ABSTRACT 


Equilibrium  configurations  for  linear  and  closed  (rings  and  regular 
polyhedra  containing  a  single  pore)  arrays  of  identical  particles  (cylinders  or 
spheres)  were  determined  by  minimizing  the  array's  surface  and  grain  boundary 
energies  (related  through  Young's  equation)  with  the  assumption  that  each  par¬ 
ticle  conserves  its  mass.  As  expected,  the  lowest  free  energy  configuration 
occurs  when  the  angle  that  defines  the  surface-grain  boundary  intersection  is 
defined  by  Young's  equation.  The  change  in  free  energy  between  the  initial  and 
equilibrium  configuration  increases  with  dihedral  angle  (i.e.,  the  equilibrium 
angle).  The  effective  sintering  'stress'  (i.e.,  driving  force)  increases  with 
the  dihedral  angle,  and  decreases  to  zero  as  the  equilibrium  configuration  is 
reached.  More  significantly,  it  is  shown  that  pores  can  shrink  to  an  equi¬ 
librium  size  if  the  number  (n)  of  coordinating  particles  is  greater  than  a 
critical  value.  The  critical  pore  coordination  number  (nc)  increases  with  the 
dihedral  angle.  Only  pores  with  n  <  nc  are  thermodynamically  unstable  during 
sintering.  It  is  also  shown  that  any  mass  transport  mechanism  will  lead  to  pore 
shrinkage  while  a  connecting  path  to  the  pore  remains  open.  Sintering  'stresses' 
increase  with  decreasing  coordination  number.  It  is  also  shown  that  the  shrink¬ 
age  strain  and  strain  rate  for  closed  arrays  increase  with  the  pore  coordination 
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number,  i.e.,  arrays  containing  more  particles  will  shrink  more  and  at  a  faster 
rate  as  they  approach  their  equilibrium  configuration.  These  results  are 
related  to  a  powder  compact  assumed  to  be  formed  by  adjoining  polyhedra.  These 
relations  suggest  that  the  resultant  force  acting  on  a  particle  will  depend  on 
the  coordination  numbers  of  polyhedra  it  shares.  These  forces  can  cause  disrup¬ 
tive  rearrangement.  The  compact's  shrinkage  strain  and  strain  rate  is  expected 
to  mimic  the  average  polyhedron,  i.e.,  is  expected  to  be  inversely  proportional 
to  the  compact's  initial  relative  density. 

1.0  INTRODUCTION 

Most  sintering  theories  have  been  based  on  ideas  developed  by  Herring,* 
who  reasoned  that  local  surface  curvature  determines  local  vacancy  concentration 
and  differential  surface  curvature  motivates  mass  transport.  Herring's  thermody¬ 
namic  concept  has  been  rigorously  applied  to  models  concerning  the  different 
mass  transport  paths  associated  with  two  touching  spheres.  *  These  models  have 
focussed  on  the  kinetics  of  mass  transport,  and  they  have  been  extended,  with 
some  success,  to  estimate  the  densification  behavior  of  powder  compact s.4,5 
Conclusions  drawn  from  these  theories  are  not  always  in  agreement  with  experi¬ 
ments.  For  example,  theory  suggests  that  densification  kinetics  should  increase 
with  decreasing  particle  size.  Yet,  common  experience  indicates  that  powders 
with  a  very  small  crystallite  size  can  be  difficult  to  density.  It  is  now  com¬ 
monly  accepted  that  agglomerates,  common  to  chemically  derived  powders,  are  one 
cause  for  poo*'  sinterability,  and  that  particle  packing  has  a  great  influence  on 
densification.  Mounting  evidence  suggests6  that  the  two-particle  models,  which 
lack  pores,  are  insufficient.  It  must  be  concluded  that  since  the  practical  aim 
of  densification  is  to  rid  the  powder  compact  of  its  void  phase,  thus,  the  void 
phase  produced  by  particle  arrangements  must  be  studied  along  with  the  mass 
transport  produced  by  the  particles. 

7  ft 

Early  studies  ’  concerned  with  the  kinetics  of  mass  transport  driven 
by  differential  surface  curvature  in  a  powder  compact  assumed  that  the  void 
phase  could  be  described  during  a  latter  stage  of  densification  as  spherical 
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pores.  Pores,  however,  are  only  spherical  in  a  polycrystalline  material  if  they 
reside  within  a  single  grain  where  the  surface  energy  is  isotropic.  As  first 
pointed  out  by  Cobley  and  further  developed  by  Kingery  and  Francois,*0  isolated 
pores  intersected  by  grain  boundaries  will  exhibit  a  surface  curvature  that 
depends  on  the  number  of  intersecting  grain  boundaries  (or,  the  number  of  grains 
coordinating  the  pore)  and  the  dihedral  angle  (defined  through  Young's  equation* 
by  ^he  ratio  of  the  grain  boundary  to  surface  energy  ratio).  Using  geometrical 
arguments,  they  showed  that  if  the  dihedral  angle  was  to  be  maintained  at  each 
grain-boundary/pore  surface  intersection,  then  the  curvature  of  the  pore  surface 
between  these  intersections  depends  on  the  pore's  coordination  number.  They 
concluded  that  if  the  coordination  number  is  less  than  a  critical  value,  then 
the  pore's  surface  will  be  concave  (looking  from  within),  which  will  promote 
mass  transport  into  the  pore  and  pore  closure.  They  also  concluded  that  if  the 
pore's  coordination  number  is  greater  than  a  critical  value,  its  surface  curva¬ 
ture  will  be  convex,  promoting  mass  transport  away  from  the  pore,  and  thus,  pore 
growth.  The  critical  pore  coordination  number  depends  on  the  dihedral  angle.  A 
more  rigorous  energy  analysis1*  has  shown  that  the  isolated  pore  will  be  stable 
(neither  shrink  nor  grow)  when  its  coordination  number  is  greater  than  the 
critical  value. 

Hoge  and  Pask*2  were  the  first  to  use  thermodynamics  to  analyze  the 
void  phase  within  a  powder  compact  by  determining  the  lowest  free  energy  con¬ 
figuration  developed  by  periodic  arrangements  of  identical  spheres.  The  spheres 
were  allowed  to  penetrate  one  another  and  conserve  their  mass  (grain  boundaries 
formed  during  penetration  remain  flat).  Invoking  thermodynamics  by  using 
Young's  equation*  they  showed  that  the  array  can  reach  an  equilibrium  configu¬ 
ration  without  complete  pore  closure  if  the  dihedral  angle  for  the  material  in 
question  is  less  than  a  critical  value.  Although  Hoge  and  Pask  did  not  develop 

a  general  relation  for  pore  closure,  they  were  the  first  to  suggest  that  under 
certain  conditions,  the  void  phase  can  form  pores  with  an  equilibrium  size. 


♦See  Eq.  (2). 
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The  packing  of  particles,  such  as  identical  spheres  or  otherwise, 
rarely  results  in  a  periodic  arrangement.  The  random,  dense  packing  of  identi- 
cal  spheres  has  been  described-1  as  an  arrangement  of  different,  irregular  poly- 
hedra  where  touching  spheres  are  centered  at  vertices.  The  different  polyhedra 
are  joined  at  corners,  edges  and/or  faces  by  shared  vertices  (spheres).  The 
void  space  within  each  polyhedron,  which  may  develop  into  a  pore,  has  three 
pertinent  properties,  viz.,  a  surface  area,  a  volume  and  a  coordination  number 
(number  of  spheres  (or  vertices)  that  define  the  polyhedron).  Using  this  idea, 
one  might  approximate  any  powder  compact  with  an  arrangement  of  polyhedra,  each 
containing  a  pore.  The  compact's  macroscopic  (bulk  scale)  behavior  can  be 
approximated  by  the  behavior  of  the  average  polyhedron.  The  coordination  number 
of  the  average  polyhedron  will  be  inversely  proportional  to  the  compact's  bulk 
density.  Microscopic  (particle  scale)  behavior  can  be  approximated  by  the  dif- 
ferential  behavior  of  the  polyhedra.  With  this  concept  in  mind,  it  would  be 
instructive  to  examine  the  behavior  of  polyhedra  which  might  be  used  to  estimate 
the  behavior  of  a  powder  compact. 

The  object  here  is  to  determine  the  behavior  of  symmetric  arrays  of 
identical  cylinders  or  spheres  that  are  allowed  to  undergo  mass  transport  to 
minimize  their  free  energy  without  losing  their  individuality,  i.e.,  the  mass  of 
each  particle  (cylinder  or  sphere)  will  be  held  constant.  The  case  where  intra- 
particle  diffusion  is  allowed,  i.e.,  some  particles  grow  as  others  disappear, 
will  be  presented  elsewhere.14  Both  linear  arrays  and  closed  arrays  (arrays 
containing  a  single  'pore')  are  examined.  Calculations  will  involve  the 
particle  free  energy  as  a  function  of  configurational  change,  determination  of 
the  array's  equilibrium  configuration,  conditions  for  pore  disappearance  and 
stability,  sintering  'stresses',  and  the  array's  shrinkage  strain.  The  results 
of  these  calculations  are  consistent  with  some  of  Cannon's15  theoretical 
observations  and  lead  to  several  interesting  conclusions  concerning  sintering 
behavior  not  taught  by  current  theory. 
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2.0  ENERGY  CALCULATIONS 


2.1  General  Approach 

Arrays  where  identical  particles  (spheres  and  cylinders)  are  identi¬ 
cally  related  to  one  another,  viz.,  lines,  rings  and  regular  polyhedra  were 
considered.  The  energy  per  particle  determined  from  the  sum  of  the  energies 
associated  with  its  surface  area  (As)  and  its  grain  boundary  area  (Ab)  formed 
with  other  adjacent  particles  is  expressed  as 

E  =  YsAs  +  YbAb  »  (1) 

where  y$  and  yb  are  the  energies  per  unit  area  associated  with  the  surface  and 
the  grain  boundary,  respectively;  both  are  assumed  to  be  isotropic.  All  other 
energies  associated  with  the  particle  are  assumed  to  be  constant  as  the  array 
changes  its  configuration  during  mass  transport.  It  was  assumed  that  ys  and  yb 
are  related  through  Young's  equation 

Yb  =  2Ys  cos  Ue/2)  ,  (2) 

where  <j»e  is  the  equilibrium  (dihedral)  angle  formed  where  the  grain  boundary 
intersects  the  particle's  surface.  Equation  (1)  can  be  rewritten  as 

E  =  ys^As  +  2  cos  (<V2^Ab^  * 

By  differentiating  Eq.  (3),  it  can  be  seen  that  the  contact  area  will  only  grow 
(dAb  >  0)  when  the  external  surface  area  decreases  (dA$  <  0),  and  it  will  cease 
to  grow  when  dAs/dAb  =  -2  cos(4»e/2),  i.e.,  when  dE  =  0.  The  configuration  of 
the  array  when  the  contact  area  ceases  to  grow  is  defined  as  the  equilibrium 
configuration.  As  defined,  the  development  of  the  equilibrium  configuration 
precludes  mass  transport  between  particles,  i.e.,  particles  maintain  their 
initial  mass  during  sintering. 
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To  determine  the  energy  of  the  array  as  the  contact  area  grows  between 
the  touching  particles,  one  needs  to  define  an  initial  particle  shape,  how  the 
particle  changes  its  shape,  and  the  relation  between  the  external  surface  area 
and  the  grain  boundary  area  during  mass  transport.  Theory  and  experiments  con¬ 
cerning  two  touching  spheres  show  that  mass  is  transported  to  the  contact  posi¬ 
tion  due  to  the  negative  surface  curvature  at  this  position  relative  to  other 
positions  on  the  particle's  surface. ^  Since  mass  transport  is  a  kinetic 
phenomenon,  the  surface  curvature  will  depend  on  the  distance  from  the  grain 
boundary  developed  between  the  particles  during  the  attempt  to  develop  an 
equilibrium  configuration.  The  nonuniform  surface  curvature  would  make  the 
definition  of  the  particle  shape  change  analytically  intractable.  Therefore,  to 
allow  analytical  calculations,  it  is  assumed  that  the  surface  curvature  is 
independent  of  position  as  the  array  develops  its  equilibrium  configuration, 
e.g.,  the  initial  touching  spheres  remained  spherical  as  they  interpenetrate, 
but  increase  their  radius  to  conserve  mass.  Relative  to  the  real  system  (surface 
curvature  is  some  function  of  position),  this  assumption  will  result  in  a 
different  array  energy  (expected  to  be  lower)  during  mass  transport  towards 
equilibrium.  But,  since  the  equilibrium  configuration  is  defined  by  particle 
surface  and  grain  boundary  configurations  that  minimize  free  energy,  it  will  be 
independent  of  the  configurational  changes  prior  to  equilibrium.  Thus,  although 
the  assumed  particle  configurational  change  used  to  calculate  the  array's  energy 
will  be  different  from  the  real  system,  their  configuration  and  energy  will  be 
identical  at  equilibrium. 

By  assuming  that  the  curvature  of  the  external  surface  is  independent 
of  position  and  each  particle  conserves  its  mass,  the  relations  between  external 
and  grain  boundary  areas  and  the  determination  of  the  surface  curvature  becomes 
a  matter  of  geometry  and  trigonometry,  and  only  become  somewhat  complex  for  the 
polyhedra  where  one  must  keep  track  of  the  number  of  spheres  that  penetrate  one 
another  as  the  configuration  changes  toward  equilibrium.  The  condition  that 
curvature  is  independent  of  position  also  imposes  one  important  limitation, 
i.e.,  calculations  concerning  spheres  are  limited  to  the  cases  where  the  dis¬ 
tance  between  mass  centers  decreases.  This  limitation  is  not  imposed  upon 
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cylinders  which  will  be  used  to  investigate  the  array  energetics  for  mass 
transport  paths  where  the  distance  between  mass  centers  is  constant  (i.e.,  ex¬ 
pected  for  transport  by  surface  diffusion  and/or  evaporation/condensation). 

Several  geometric  variables  can  be  used  to  describe  the  configurational 
change  of  the  array  as  the  contact  area  grows.  They  include  1)  the  angle  de¬ 
fined  by  the  tangents  to  the  external  surface  at  the  grain  boundary,  termed  the 
contact  angle  (4.,  see  Fig.  lb),  2)  the  grain  boundary  length  (h),  and  3)  the 
radius  (Rp,  see  Fig.  4)  of  the  circle  (or  sphere)  that  circumscribes  the  grain 
boundary  intersections  on  the  internal  surfaces  and  describes  the  pore  size  for 
closed  arrays.  Although  these  variables  are  interrelated,  and  for  the  appro¬ 
priate  cases  can  be  used  interchangeably,  it  is  most  convenient  to  express  the 
configurational  changes  in  terms  of  the  contact  angle,  which  is  equal  to  zero 
(4>  =  0)  when  the  particles  are  just  touching,  and  increases  as  the  array 
develops  its  equilibrium  configuration.* 


2.2  Infinite,  Linear  Arrays 


2.2.1  Cylinders;  Approaching  Centers 


Consider  an  infinite  line  of  identical  touching  cylinders,  with  an 
initial  radius  r^ ,  shown  in  Fig.  1  at  different  stages  of  interpenetration.  As 
mass  centers  approach  and  the  cylinders  interpenetrate,  their  radius  (r)  in¬ 
creases.  Using  the  geometrical  variables  shown  in  Fig.  lb  and  Eq.  (2),  it  can 
be  shown  that  the  energy  per  unit  length  of  cylinder  can  be  written  as  a  func¬ 
tion  of  four  variables,  viz.,  the  contact  angle  (<!>),  the  normalized  radius 
(R  =  r/r.j),  normalized  grain  boundary  length  (H  =  h/r^),  and  the  dihedral  angle 
( 4»e ) : 


E 


2wV< 


R  ( n  - 


4-) 


+  H  cos  Ue/2) 


1 

u 


(4) 


*As  shown  in  Appendix  1,  the  equilibrium  configuration  can  also  be  determined 
by  a  Wulff  construction. 
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Since  each  particle  maintains  its  initial  mass,  one  can  define  H  in  terms  of 
R  and  4»: 


u  _  tc/R  +  R ( 41  -  't  +  sin  4>) 

H  l  cos  4>/Z  * 

Substituting  Eq.  (5)  into  (4),  one  obtains 
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(6) 


The  equilibrium  condition  is  determined  by  minimizing  particle  energy  with 
respect  to  the  geometrical  variables  4),  and  R.  It  can  be  shown  that  when 
dE(<|»,R)  =  0,  H  =  2r  sin  4>/2,  which  allows  us  to  redefine  Eq.  (6)  as 


and 


E  =  2uVs 


R  =  (- 


n  -  4>  +  2  sin  y  cos  — 

y  n(n  -  +  sin  <|j) 

it  'i  1/2 


u  -  4<  +  sin  <|> 


(7) 

(a) 


Figure  2  shows  that  the  particle's  energy  decreases  with  contact  angle  and 
reaches  a  minimum  at  (|<  =  4>e.  Also,  lower  energies  result  for  larger  dihedral 
angles. 

Particle  energy  can  also  be  related  to  the  distance  of  particle- 
particle  approach,  which  can  be  expressed  in  units  of  strain,  i.e.,  e  =  1  -  R 
cos  4^/2 .  As  shown  in  Fig.  3  particle  energy  decreases  with  particle-particle 
approach  and  reaches  a  minimum  energy  condition  at  increasingly  larger  strains 
with  increased  dihedral  angle  (4>e).  (Shrinkage  strains  will  be  discussed  fur¬ 
ther  in  Section  4.) 
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2.2.2  Cylinders;  Center-to-Center  Distance  Constant 

When  the  particle  centers  of  mass  are  fixed  the  radius  of  curvature  of 
the  external  surface  (r)  is  easily  determined  from  geometry  (i.e.,  r  =  ^/cosU/2) 
The  centers  of  curvature  are  neither  coincident  with  one  another  nor  with  the 
center  of  mass. 

It  can  be  shown  that  the  energy  per  unit  length  can  be  written  as  a 
function  of  the  contact  angle  (4>)  and  the  dihedral  angle  ( <Pe ) : 


E  =  2ur.  y 
i  Ts 


(1  *  -  4*  -  Sin  4n  e  x 

i  - 21  > cos  ir 


2n  cos 


n  -  iji 

d / 

7i  cos  y 


(9) 


The  equilibrium  configuration  occurs  when  dE/d4>  =  0,  which  can  be  shown  to  occur 
when  4>  =  4>g. 

Figure  4  shows  particle  energy  for  both  approach  (solid  lines,  Eq.  (7)) 
and  nonapproach  (broken  lines,  Eq.  (9))  of  particle  centers.  As  shown,  energies 
are  lower  for  the  case  of  particle-particle  approach. 

2.2.3  Spheres 

For  the  case  of  spheres,  the  restriction  that  surface  curvature  is 
independent  of  position  limits  the  calculation  to  conditions  were  the  mass 
centers  approach  one  another.  (The  case  where  the  distance  between  mass  centers 
is  fixed  results  in  barrel -shaped  particles.)  Particle  radius  of  curvature  (r) 
is  therefore  uniquely  defined  by  the  particle's  initial  radius  (r^ )  and  the 
contact  angle  (4>): 


r  =  n  ^  cos  y  (2  +  sin2  j  ^ 
The  energy  per  particle  can  be  expressed  as: 


(10) 


4«r(2Ys 


4.  1  .24,  ^e 

cos  y  *  Tg  sin  y  cos  — 

4  (2  +  sin2 


^  cos  ^  ^  t  iUI 
and  is  a  function  of  the  contact  angle  (<|>),  and  dihedral  angle  ( 4>e ) . 
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Similar  to  the  case  of  the  linear  array  of  cylinders,  the  energy  of  a 
particle  decreases  and  reaches  equilibrium  when  4>  =  4<e.  The  change  in  energy 
for  different  values  of  4>e  is  similar  to  that  shown  in  Fig.  2  for  the  cylinders. 


2.3  Closed  Arrays 

Forming  a  ring  array  from  a  linear  array,  or  a  polyhedron  from  a  planar 
array,  introduces  several  new  geometrical  variables  needed  for  array  calcula¬ 
tions.  These  new  variables  include  the  radius  ( Rp )  of  the  pore  enclosed  by  the 
particles,  and  the  finite  number  (n)  of  particles  that  coordinate  the  pore.  The 
pore  radius  can  be  defined  by  a  circle  (or  sphere)  that  circumscribes  the  posi¬ 
tions  where  the  grain  boundaries  intersect  the  pore's  surface  as  shown  in 
Fig.  5.  The  radius  of  the  pore  is  dependent  on  other  array  variables.  The 
pore's  coordination  number  is  an  independent  variable  and  can  be  described  by 
the  angle  that  subtends  a  single  particle.  For  ring  arrays  the  coordination 
angle  is  simply  related  to  the  coordination  number,  0  =  2n/n. 

One  other  important  difference  is  that  the  center  of  mass  for  particles 
in  closed  arrays  moves  during  mass  transport  (e.g.,  particle  cross  sections 
become  'pie'  shaped),  and  the  center  of  curvature  is  no  longer  coincident  with 
the  center  of  mass.  This  observation  is  more  important  in  a  latter  section 
concerning  shrinkage  strains. 

2.3.1  Ring  of  Cylinders;  Approaching  Centers 

It  can  be  shown  that  when  <\>  +  2n/n  <  u,  the  energy  per  particle  in  the 
ring  array  is  identical  to  that  for  the  linear  array,  viz.,  as  expressed  by 
Eq.  (7).  For  the  condition  <\>  +  2it/n  >  u  the  particle  energy  can  be  shown  to  be 
a  function  of  contact  angle  (41),  dihedral  angle  (4>e),  and  pore  coordination  (n): 

rc  +  —  -  (j,  +  (sin  i  +  cos  -y/tan  — )  2  cos  iX 

E  =  2nr .  y  - \ 1 - "1 -  *1  (12) 

1  s  2^(cos*  J/tan  £  +  £  sin  4-  + 


The  pore  size  is  given  by 


10 

J  7623TC / j  bs 


Rockwell  International 

Science  Center 


R 
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cos  $  +  j)  I  -  — 

sin  (2)  \  11  "  *  +  sln  * 


ri 


(13) 


To  determine  conditions  of  pore  stability/disappearance,  the  energy  per 
particle  (or  array  energy)  can  be  related  to  the  pore  size.  As  shown  in  Fig.  6 
(for  the  condition  4>e  =  150°)  a  minimum  energy  can  occur  before  the  pore  disap¬ 
pears,  i.e.,  further  pore  closure  will  increase  the  free  energy  of  the  array. 

The  conditions  required  for  pore  closure  can  best  be  determined  by 
substituting  <]>  -  4>e  in  Eq.  (13)  to  obtain  the  relation  for  the  equilibrium  pore 
size.  Examination  of  Eq.  (13)  shows  that  the  pore  will  have  a  finite  size  when 
cl>e  +  2u/n  <  u  and  will  disappear  when  +  2n/n  >  n.  The  condition  where 
4<e  +  2n/n  =  n  defines  the  critical  coordination  number,  viz.,  nc  =  2n/(it  -  4*e). 
Pores  will  disappear  when  n  <  nc  and  shrink  to  an  equilibrium  size  when  n  >  nc. 

2.3.2  Ring  of  Cylirders;  No  Mass  Transport  Between  Pore  and  External  Surface 

With  mass  transport  prevented  between  the  inside  pore  surface  and  the 
exterior  surface  of  the  array,  material  which  is  outside  the  polygon  formed  by 
connecting  the  centers  of  adjoining  particles  (Fig.  5a)  cannot  be  transported  to 
the  pore  surface.  This  requires  that  the  inside  and  outside  regions  separately 
conserve  mass.  The  equilibrium  configuration  given  these  mass  transport  con¬ 
straints  will  be  determined,  via  Eq.  (3),  by  calculating  A$  and  Ab  with  the 
assumption  that  the  interior  and  the  exterior  surfaces  of  the  array  have  uniform 
curvatures. 

It  can  be  shown  that  the  radius  of  curvature  of  the  interior-pore  (R^  ) 
and  exterior-surface  (R0)  are  a  function  of  the  pore  coordination  (n)  contact 
angle  (<|>)  and  initial  particle  radius  ( r^ )  and  generally  are  not  equal: 

2n  0  .  n 

w  -  - - 2  cot  — 

_ n _ n _ _ 

it - 5  -  <j,  +  sin  <j>  -  2  cos  $  cot  — 

n  c  n 


X' 
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_ _ +A 


\Km 
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where  is  (  +  )  for  +  /p  <  u 


R  =  r.K 
o  i  o 


K  = 
o 


u  +  ~  +  2  cot  £ 
n  n 


a/2 


(14b) 


n  +  —  -  4-  +  sin  4>  +  2  cos2  4"  cot  — 
n  c  n 


Particle  energy  can  also  be  decoupled  into  interior  and  exterior  parts, 
and  these  parts  have  been  arbitrarily  separated  by  the  polygon  formed  by  con¬ 
necting  particle  centers.  The  energy  of  the  interior  pore  region  (E^)  was  taken 
to  consist  of  both  the  pore  surface  energy,  and  the  grain  boundary  energy  lying 
within  the  polygon  formed  by  connecting  particle  centers,  while  the  energy  of 
the  exterior  region  (E0)  was  taken  to  consist  of  both  the  external  surface 
energy  and  the  grain  boundary  energy  lying  outside  the  polygon.  The  energy  of 
the  interior  (E ^ )  and  exterior  (EQ)  region,  per  unit  length,  is  a  function  of 
the  pore  coordination  number  (n),  contact  angle  (<i>),  and  dihedral  angle  ( 4>e ) : 


E.  = 


r.K. 
l  1 


4 •  ~  zr  ~  2  cos  +  ■£) 


+  2r . 


<1> 

cos 

tan  — 
n 


(15a; 


E  = 


o  ,  cos 

£+2c°5 


4) 
e 

cos  "2 — 


ji 

n 


cos 


-  2ri 


tan  — 
n 


(15b) 


The  equilibrium  configuration,  determined  with  d(E^  +  E0 ) /d 4>  =  0, 
occurs  when  4»  =  «|»e.  The  curvatures  of  the  interior  pore  surface  and  exterior 
surface  for  the  equilibrium  configuration  can  be  determined  by  substituting  4>e 
for  4>  in  Eqs.  (14a)  and  (14b).  The  two  curvatures  are  only  equal  when  their 
centers  of  curvature  are  coincident.  For  4>e  +  2u/n  <  u  the  radius  of  curvature 
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of  the  interior  pore  surface  is  positive  (convex  as  viewed  from  the  pore).  For 
<j>e  +  2u/n  >  n,  R^  is  negative  (concave).  The  radius  of  curvature  of  the  ex¬ 
ternal  surface  (RQ)  is  always  positive  (convex  as  viewed  from  outside  the  ring). 
In  all  cases,  the  radius  of  curvature  of  the  internal  surface  is  always  such, 
that  if  mass  could  be  transported  between  the  two  surfaces,  the  pore  would 
either  shrink  or  disappear  as  described  in  the  previous  section  (mass  centers 
approach) . 

2.3.3  Ring  of  Spheres 

Rings  of  spheres  exhibit  similiar  behavior  as  rings  of  cylinders.  The 
requirement  of  uniform  surface  curvature,  however,  restricts  the  calculations  to 
particle  approach.  For  conditions  where  4>e  +  2u/n  <  it,  it  can  be  shown  that  the 
particle  energy  is  identical  to  that  determined  for  the  linear  array  (Eq.  11); 
particle  energy  for  these  conditions  is  not  a  function  of  pore  coordination 
number.  For  conditions  where  the  pore  disappears  (ii>e  +  2u/n  >  n),  all  particles 
penetrate  one  another  and  the  determination  of  the  particle  energy  is  somewhat 
complex  and  given  in  Appendix  2. 

The  conclusions  concerning  the  ring  of  spheres  are  similar  to  the  ring 
of  cylinders  with  approaching  centers,  i.e.,  when  4>e  +  2n/n  <  n,  the  equilibrium 
configuration  contains  a  stable  pore,  and  when  <\>Q  +  2n/n  >  the  pore 
disappears. 

2.3.4  Regular  Polyhedra  (See  Appendix  2  for  Details) 

Only  polyhedra  consisting  of  identical  spheres  symmetrically  arranged 
can  form  equilibrium  structures  consisting  of  surfaces  of  uniform  curvature. 

This  limits  our  analysis  to  the  five  regular  polyhedron,  as  described  in  Table  2 
of  Appendix  II.  The  energy  calculations,  more  complex  due  to  the  varying 
geometry,  are  detailed  in  Appendix  II. 

It  can  be  shown  that  conclusions  concerning  equilibrium  configurations 
and  pore  stabi 1 ity /disappearance  are  similar  to  that  of  the  ring  configurations 
where  centers  approach. 
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3.0  SINTERING  'STRESSES' 


It  is  often  stated  that  the  change  in  surface  energy  is  the  driving 
force  for  sintering.  More  specifically,  the  differential  of  surface  energy  with 
respect  to  either  displacement  or  volume  dimensionally  corresponds  to  either 
force  or  stress,  respectively.  Although  these  'forces'  and  'stresses'  have  no 
mechanical  origin,  they  can  manifest  themselves  in  mechanical  action.  For 
example,  when  a  powder  compact  densities  around  an  inclusion,  the  sintering 
stresses  are  manifested  as  mechanical  stresses  which  arise  within  the  inclusion 
and  the  surrounding  powder  compact  that  can  cause  crack  extension  within  the 
powder  compact.17  From  a  densification  standpoint,  the  magnitude  of  the 
sintering  'stress'  should  defined  the  system's  potential  for  densification. 

In  the  above  section,  the  energy  of  different  arrays  was  determined  as 
a  function  of  their  geometric  changes  produced  by  intraparticle  mass  transport. 
Using  these  results  it  is  a  simple  matter  to  redefine  the  energy  in  terms  of  the 
distance  between  mass  centers  and  then  differentiate  to  obtain  the  sintering 
'force'  acting  between  the  mass  centers.  The  sintering  'stress'  can  then  be 
determined  by  dividing  the  force  by  the  grain  boundary  area  over  which  the 
normal  force  is  acting.  For  closed  particle  arrays  one  can  determine  the 
apparent  hydrostatic  pressure  which  produces  an  equivalent  sintering  force. 

3.1  Linear  Arrays 

3.1.1  Cyl i nders 

As  previously  determined  in  Section  2.2.1,  particle  energy  of  a  linear 
array  of  cylinders  is  expressed  as  Eq.  (7).  The  distance  between  mass  centers 
is  given  by 

1/2 


x  =  cos  U/2)  sW  „,) 

The  sintering  'force,'  calculated  as 


(16) 
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F  =  (dE/d<|>)/(dx/d<j>)  , 


is  given  in  the  following  expression  for  the  linear  array  of  cylinders 

F  =  2ys(cos  <|»/2  -  cos  4>e/2 )  J^-2-^  +  cot  q»/2j  .  (18) 

The  sintering  ‘stress',  per  unit  length  of  cylinder,  determined  by  dividing 
Eq.  (18)  by  the  grain  boundary  area  per  unit  length  of  cylinder  [2r  sin(^/2)], 
is  given  by 


,  (b  ^e  \  f  2  ,  <j>\  1 

*c  =  77  (cos  7  '  cos  cot  7 i  “ — $ 


Tt~  4>  +  sin  (|i  ( 19 ) 


As  shown  by  Eq.  (19),  the  sintering  'stress'  is  inversely  proportional  to  the 
initial  particle  size  (r^).  As  shown  in  Fig.  7,  it  decreases  to  zero  as  4,  ■*  4>e, 
and  increases  with  increasing  dihedral  angle. 

3.1.8  Spheres 

In  a  similar  manner,  it  can  be  shown  that  the  sintering  'stresses'  for 
the  linear  array  of  spheres,  can  be  expressed  as 


Yc  .  4-  1  +  cos2  $ 


Js  "  TT 


(cos  y  -  COS  y-) 


2  '  !  8  4>  ijj  ’ 

sin  y  cos  f 


and  exhibits  the  same  behavior  previously  described  for  a  linear  array  of 
cylinders. 

3.2  Closed  Arrays 

The  apparent  hydrostatic  sintering  'pressure'  was  only  determined  for 
the  ring  array  of  cylinders.  With  knowledge  of  the  'force'  acting  normal  to  the 
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grain  boundaries  (Eq.  (18)),  one  can  determine  an  apparent  hydrostatic  pressure 
acting  on  the  array  in  the  same  way  as  one  determines  the  stress  in  a  thin  wall 
cylinder.  That  is,  the  apparent  hydrostatic  pressure  is  determined  by  equating 
the  component  of  the  force  (F),  in  the  radial  direction,  to  the  hydrostatic 
pressure,  multiplied  by  the  component  of  surface  area  perpendicular  to  the 
radial  direction.  The  apparent  hydrostatic  pressure  (per  unit  length  of 
cylinder)  is 


P 


sin 


n 

n 


cos 


V 


u  -  +  sin  (\> 

n 


(21) 


Equation  (21)  shows  that  P  is  inversely  related  to  the  initial  particle  size 
(r.j)  and  decreases  to  zero  at  the  equilibrium  condition  (<|>  =  ^e).  Also,  as 
shown  in  Fig.  8,  P  increases  with  decreasing  coordination  number. 

Although  sintering  'pressures'  were  not  detailed  for  the  polyhedra,  an 
examination  of  their  energy  relations  shows  that  similar  results  would  be 
obtained  as  shown  in  Fig.  8. 


4.0  SHRINKAGE  STRAINS 

Engineering  strains  are  usually  used  to  define  the  shrinkage  behavior 
of  a  powder  compact,  thus  shrinkage  strains  determined  here  are  related  to  the 
initial  array  dimension. 

Array  shrinkage  must  be  related  to  the  particle's  center  of  mass,  which 
for  a  linear  array,  is  identical  to  the  center  of  curvature.  For  closed  arrays 
the  center  of  mass  changes  its  position  with  respect  to  the  center  of  curvature 
as  the  array  changes  its  configuration  (see  Fig.  5). 
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4.1  Linear  Array 

The  shrinkage  strain  for  linear  arrays  of  cylinders  and  spheres  is  only 
a  function  of  contact  angle  ( 4* )  as  shown  by 


ecyl 


_ U/2 

4/  +  sin  4»' 


(22a) 


'sph 
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2  cos2  |  1/3 


2  +  sin 


11 
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(22b) 


The  maximum  shrinkage  occurs  when  the  array  reaches  its  equilibrium  configura¬ 
tion  (<|>  =  4>  ).  The  equilibrium  shrinkage  strain  increases  with  dihedral  angle 

ue). 


4.2  Closed  Arrays 

Center  of  mass  calculations  of  particle  arrays  containing  spheres  are 
more  tedious  than  for  cylindrical  arrays,  and  do  not  lead  to  additional  physical 
insight.  For  this  reason,  only  rings  of  cylinders  were  treated  in  detail. 

It  can  be  shown  that  prior  to  pore  closure  (i^+  2u/n  <  u)  the  shrinkage 
strain  for  the  ring  of  cylinders  is  given  by 


=  1  -  (cos 


4> 
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sin  — 
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4>  +  sin  4> 
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and  after  pore  closure  (4<  +  2u/n  >  n)  by 
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u  -  4>  +  —  +  sin  +  2 


tan  — 


As  shown  in  Fig.  9,  shrinkage  strain  increases  with  the  pore  coordina¬ 
tion  number.  Strains  incorrectly  calculated  from  the  center  of  curvature,  on 
the  other  hand,  do  not  show  any  variation  with  pore  coordination  number. 
Negligible  but  finite  array  shrinkage  occurs  after  pore  closure,  due  to  the 
further  cyl i nderi zation  of  the  ring  array. 


5.0  DISCUSSION 

The  thermodynamic  approach  to  the  sintering  of  simple  particle  arrays 

has  both  quantified  some  expected  behavior  and  produced  some  unexpected  results. 

It  was  shown  that  the  energy  change  and  the  driving  force  for  densifi- 
cation  increased  with  the  dihedral  angle.  This  result  is  expected,  and  suggests 
that  sintering  aids  that  decrease  the  grain  boundary  energy,  and  thus  increase 
the  dihedral  angle,  will  increase  the  driving  force  for  sintering. 

More  important,  the  analysis  for  closed  arrays  showed  that  pore  closure 

is  not  merely  a  matter  of  kinetics  as  usually  postulated.  Sound  thermodynamic 
reasoning  shows  that  pores  can  shrink  to  an  equilibrium  size.  Identical  results 
are  obtained  for  isolated  pores  in  large,  but  finite  polycrystalline  bodies.11 
Both  analyses  show  that  pore  closure  depends  on  the  pore's  coordination  number 
and  the  dihedral  angle.  Pores  shrink  to  an  equilibrium  size  when  the  sum  of  the 
dihedral  angle  and  the  coordination  angle  is  <  n  (n  >  nc)  and  disappear  when 
this  sum  is  >  it  (n  <  nc).  Similar  results  were  obtained  for  a  planar  array  by 
Cannon1^  who  used  a  different  approach  to  define  equilibrium  configurations. 

Since  the  critical  coordination  number  nc  increases  with  the  dihedral 
angle,  sintering  aids  that  increase  the  dihedral  angle  will  allow  more  pores 
within  a  given  powder  compact  to  be  thermodynamically  unstable  and  disappear. 

The  concept  that  some  pores  can  be  thermodyamically  stable  leads  one  to  hypothe¬ 
sis  ways  of  reducing  their  coordination  number.  As  detailed  elsewhere, ^  grain 
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growth  will  reduce  the  pore's  coordination  number  and  cause  stable  pores  to 
become  unstable  and  disappear.  Thus,  grain  growth,  once  thought  to  be  the 
scourge  of  sintering  phenomena,  can  be  helpful.  The  thermodynamics  of  pore 
closure  also  describes  the  role  of  particle  packing.  Namely,  compacts  with 
higher  bulk  densities  are  expected  to  contain  fewer  highly  coordinated  pores, 
and  thus,  for  a  given  sintering  schedule,  a  higher  end-point  density. 

It  should  be  noted  that  one  did  not  need  to  define  the  mass  transport 
path  for  calculations  concerning  closed  arrays  with  approaching  centers.  As 
long  as  there  remains  an  open  path  between  the  external  surface  of  the  array  and 
the  pore  surface,  surface  diffusion  and  evaporation/condensation  will  lead  to 
the  same  configurational  changes  that  are  usually  only  associated  with  transport 
mechanisms  that  cause  mass  centers  to  approach  for  linear  arrays.  That  is, 
transport  mechanisms  that  cannot  lead  to  the  shrinkage  of  linear  arrays  can  re¬ 
sult  in  the  shrinkage  of  closed  arrays.  It  must  be  concluded  that  although  sur¬ 
face  diffusion  and/or  evaporation/condensation  will  not  produce  bulk  shrinkage 
in  a  powder  compact,  some  pores  (those  with  a  lower  than  average  coordination 
number)  will  shrink  and  possibly  disappear,  while  others  (those  with  a  greater 
than  average  coordination  number)  will  grow,  i.e.,  these  transport  mechanisms 
can  lead  to  local  densification. 

It  was  shown  that  the  shrinkage  strain  for  the  closed  arrays  increases 
with  coordination  number.  This  result  is  not  unexpected  since  it  is  common  ex¬ 
perience  that  the  shrinkage  of  powder  compacts  is  inversely  proportional  to 
their  initial  bulk  density.  The  following  argument  will  show  that  the  shrinkage 
strain  rate,  which  cannot  be  explicitly  calculated  without  a  kinetic  model,  will 
also  increase  with  the  pore  coordination  number.  This  argument  is  based  on  the 
fact  that  as  the  conf iguration  changes  to  the  point  where  the  pore  disappears, 
the  mass  transport  path  is  independent  of  the  coordination  number,  i.e.,  the 
transport  path  is  identical  from  particle  to  particle  and  array  to  array.  That 
is,  within  the  same  period,  the  amount  of  shrinkage  will  depend  on  the  array's 
coordination  number,  but  the  amount  of  mass  transport  per  particle  will  be  the 
same  for  all  arrays.  Since  within  the  same  period,  the  amount  of  shrinkage  is 
proportional  to  the  arrays'  coordination  number,  both  the  shrinkage  strain  and 
strain  rate  of  an  array  will  be  proportional  to  its  coordination  number.  Thus, 


19 

07623TC / jbs 


Rockwell  International 

Science  Center 


an  array  with  a  larger  coordination  number  will  shrink  more  and  at  a  greater 
rate  than  an  array  with  a  smaller  coordination  number. 

Using  the  shrinkage  behavior  of  closed  arrays,  we  can  now  estimate  the 
behavior  of  powder  compacts.  As  suggested  above,  the  coordination  number  of  the 
average  polyhedron  that  describes  the  bulk  behavior  will  be  inversely  propor¬ 
tional  to  the  compact's  bulk  density.  One  might  expect  that  as  pores  shrink  to 
their  equilibrium  size,  the  compact's  shrinkage  strain  and  strain  rate  will  be 
inversely  proportional  to  its  initial  bulk  density.  This  behavior  has  been 
observed  in  sintering  studies  where  the  same  powder  was  compacted  to  different 
initial  bulk  density.1*1 

As  indicated  above  the  shrinkage  strain  and  strain  rate  will  depend  on 
the  pore's  coordination  number.  Since  different  polyhedra  will  shrink  differ¬ 
ently,  one  could  expect  the  differential  behavior  of  connected  polyhedra  will 
give  rise  to  stresses  during  sintering. 

Rearrangement  is  the  term  commonly  used  to  describe  the  nonuniform 
movement  of  particles  during  the  initial  stages  of  sintering.  Exner's  sintering 
studies16  of  random,  planar  arrays  of  identical  spheres  graphically  illustrate 
this  phenomenon.  Some  groups  of  particles  shrink  upon  themselves,  forming 
large,  irregular  voids  between  the  shrinking  groups.  Mercury  infiltration 
experiments  on  powder  compacts6,19  also  reveal  this  phenomena.  Small  equivalent 
capillaries  close  during  the  early  stages  of  sintering,  and  concurrently,  larger 
capillaries  grow  even  larger.  As  detailed  above,  the  magnitude  of  the  radial 
component  of  the  sintering  'force'  acting  on  a  particle  within  a  close  array  is 
inversely  proportional  to  the  coordination  number.  Particles  that  share  two  or 
more  polyhedra  within  a  powder  compact  therefore  can  have  a  resultant  force  that 
pulls  it  into  the  polyhedron  with  the  smallest  coordination  number,  whereas  a 
neighboring  particle  could  be  pulled  in  the  opposite  direction.  On  the  particle 
scale,  it  can  be  seen  that  forces  acting  on  one  particle  will  be  different  from 
particle  to  particle  and  will  certainly  give  rise  to  rearrangement.  On  the 
polyhedra  scale,  polyhedra  with  coordination  numbers  greater  than  the  average 
will  try  to  shrink  more  and  at  a  greater  rate  than  the  average.  Constraint  by 
the  average  shrinkage  behavior  will  give  rise  to  disruptive,  tensile  stress.  On 
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the  agglomerate  scale,  it  has  been  experimentally  shown^  that  agglomerates  with 
a  lower  bulk  density  than  the  average  will  produce  circumferential,  crack-like 
voids.  Thus,  disruptive  processes  due  to  differential  shrinkage  behavior  are 
expected  to  occur  at  all  scales.  Voids  produced  by  disruptive  processes  will 
have  higher  coordination  numbers  than  those  present  before  disruption.  Thus, 
rearrangement  processes  due  to  nonuniform  packing  can  be  the  most  detrimental 
phenomenon  that  occurs  during  sintering. 

Periodic  packing,  which  requires  identical  particles,  might  be  sug¬ 
gested  as  a  packing  arrangement  that  would  eliminate  resultant  particle  forces, 
and  thus  disruptive,  rearrangement  phenomena.  Unfortunately,  disruptive,  resul¬ 
tant  forces  would  still  exist  for  particles  located  at  boundaries  between  order 
domains,  which  could  produce  crack-like  voids  as  large  as  the  domain  size. 

Thus,  unless  periodic  packing  could  be  obtained  without  domain  boundaries,  one 
might  conclude  that  packing  arrangements  such  as  those  based  on  modifi cations  of 
dense  random  packing  might  minimize  the  size  of  the  disruptive  unit. 


ACKNOWLEDGEMENTS 

This  work  was  supported  by  the  Office  of  Naval  Research  under  Contracts 
No.  N00014-84-K-0286  at  IJCLA,  and  No.  N00014-84-C-0298  at  Rockwell  International 
Science  Center. 


6.  REFERENCES 

1.  C.  Herring,  "Surface  Tension  as  a  Motivation  for  Sintering,"  The  Physics 
of  Powder  Metallurgy,  ed.  T.E.  Kii  gston,  Chapt.  8,  McGraw-Hill,  NY  (1951). 

2.  W.D.  Kingery  and  M.  Berg,  "Study  of  the  Initial  Stages  of  Sintering  Solids 
by  Viscous  Flow,  Evaporation-Condensation,  and  Self-Diffusion,"  J.A.P.  2b , 
1205  (1955). 


21 

J  76  23TC / j  bs 


Rockwell  International 

Science  Center 


3.  F.A.  Nichols,  "Coalescence  of  Two  Spheres  by  Surface  Diffusion,"  J.A.P. 

2]_,  2805  (1966). 

4.  J.6.R.  Rockland,  "The  Determination  of  the  Mechanism  of  Sintering,"  Acta. 
Metall .  _15,  277,  (1967). 

5.  D.  Lynn  Johnson  and  Ivan  B.  Cutler,  “Diffusion  Sintering:  I.  Initial  Stage 
Sintering  Models  and  Their  Application  to  Shrinkage  of  Powder  Compacts," 

J.  Am.  Ceram.  Soc.  46_,  541  (1963). 

6.  F.F.  Lange,  "Si nterabi 1 ity  of  Agglomerated  Powders,"  J.  Am.  Ceram.  Soc. 

67_,  83,  (1984). 

7.  R.L.  Coble,  "Sintering  Crystalline  Solids.  I.  Intermediate  and  Final  State 
Diffusion  Models,"  J.A.P.  32,  787,  (1961). 

8.  G.C.  Kuczynski,  "The  Mechanism  of  Densif ication  During  Sintering  of 
Metallic  Particles,"  Acta.  Metall.  4_,  58,  (1956). 

9.  R.L.  Coble,  "Diffusion  Sintering  in  the  Solid  State,"  in  Kinetics  of  High 
Temperature  Processes,  ed.  W.D.  Kingery,  147-163,  Joint  Publication  of  the 
Technology  Press  of  M.I.T.  and  John  Wiley  &  Sons,  Inc.,  New  York,  (1959). 

10.  W.D.  Kingery  and  B.  Francois,  "Sintering  of  Crystalline  Oxides,  I: 
Interaction  Between  Grain  Boundaries  and  Pores,"  in  Sintering  and  Related 
Phenomena,  ed.  G.C.  Kuczynski,  N.A,  Hooten  and  G.F.  Gilbon,  471-98,  Gordon 
Breach,  (1967). 

11.  B.J.  Kellett  and  F.F.  Lange,  "Thermodynamics  of  Densificiation,  Part  III: 
Experimental  Relation  Between  Grain  Growth  and  Pore  Closure,"  submitted  to 
J.  Am.  Ceram.  Soc. 

12.  C.E.  Hoge  and  J.A.  Pask,  Ceramurgia  International  3  [3],  95-9,  (1977). 


22 

J7623TC / j  bs 


Rockwell  International 

Science  Center 


H.J.  Frost,  "Overview  17:  Cavities  in  Dense  Random  Packing,"  Acta.  Metall. 
30  [5],  889-904,  (1982). 

F.F.  Lange  and  BJ.  Kellett,  "Thermodynamics  of  Densification,  Part  II: 
Grain  Growth  in  Porous  Compacts  and  Relation  to  Densification,"  submitted 
to  J.  Am.  Ceram.  Soc. 

R.M.  Cannon,  "On  the  Effects  of  Dihedral  Angle  and  Pressure  on  the  Driving 
Force  for  Pore  Growth  or  Shrinkage,"  unpublished  work,  (1981). 

H.E.  Exner,  "Principles  of  Single  Phase  Sintering,"  Rev.  Powder  Metall. 
Phys.  Ceramics  _1_,  1-251,  (1979). 

F. F.  Lange  and  M.  Metcalf,  "Processing-Related  Fracture  Origins:  II. 
Agglomerate  Motion  and  Cracklike  Internal  Surfaces  Caused  by  Differential 
Sintering,"  J.  Am.  Ceram.  Soc.  _66_  [63,  398-406,  (1983). 

B.  Kellett  and  F.F.  Lange,  "Stresses  Induced  by  Differential  Sintering  in 
Powder  Compacts,"  J.  Am.  Ceram.  Soc.  _67_,  369,  (1984). 

O.J.  Whitermore,  Jr.,  and  J.J.  Sipe,  "Pore  Growth  During  the  Initial 
Stages  of  Sintering  Ceramics,"  Powder  Technol .  9^,  159,  (1974). 

G.  Wulff,  Z.  Kristallogr.  34,  449,  (1901). 

W.L.  Winterbottom,  "Equilibrium  Shape  of  a  Small  Particle  in  Contact  With 
a  Foreign  Substrate,"  Acta.  Metall.  _15_,  303,  (1966). 

C. A.  Shell  and  N.H.  Fletcher,  "Decoration  Criteria  for  Surface  Steps," 
Acta.  Metall.  18,  1083,  (1970). 


23 

J7623TC/jbs 


Rockwell  International 

Science  Center 


* 


SC85-30277 


a)  Initial,  b)  generalized,  and  c)  equilibrium  conf iguration  of  a 
linear  array  of  cylinders,  rj  is  the  initial  particle  radius,  r  par¬ 
ticle  radius,  is  the  contact  angle  defined  by  the  tangents  of 
surface  at  the  grain  boundary,  h  is  the  grain  boundary  length  and 
is  the  dihedral  angle  defined  by  Young's  equation  (Eq.  (2)). 
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Particle  energy  per  unit  cylinder  length,  normalized  by  the  initial 
particle  energy,  for  the  linear  array  of  cylinders,  as  a  function  of 
the  contact  angle  ( u<) ,  for  different  dihedral  angles  (ne). 
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Particle  energy  per  unit  cylinder  length,  normalized  by  the  cylin¬ 
der's  initial  energy,  for  a  linear  array  of  cylinders  for  the  con¬ 
dition  of  particle-particle  approach  (dashed  line)  and  no  particle 
particle  approach  (solid  line),  as  a  function  of  the  contact  angle 
M,  for  different  dihedral  angles  (ui  ). 
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Particle  energy  per  unit  length  of  cylinder,  normalized  by  the  ini¬ 
tial  energy  of  a  cylinder,  as  a  function  of  circumscribed  pore  radius 
(R  ),  for  different  pore  coordination  numbers  (n).  Results  are 
illustrated  for  a  dihedral  angle,  ne  =  150  . 
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APPENDIX  1 

on 

As  was  shown  by  Wulff,  the  equilibrium  configuration  of  an  isolated 
particle  can  be  determined  from  knowledge  of  its  surface  energies.  Basically, 
the  Wulff  construction  determines  a  particle's  equilibrium  shape  as  the  interior 
envelope  produced  by  perpendicular  planes  placed  at  the  ends  of  the  particle's 
surface  energy  vectors.  These  vectors  emanate  from  a  single  point  (Wulff 
point),  and  have  an  orientation  perpendicular  to  the  crystallographic  planes 
they  represent.  The  Wulff  construction  of  a  particle  of  isotropic  surface 
energy  is  a  sphere. 

The  Wulff  construction  of  an  array  of  identical  particles  (or  any  array 
which  forms  flat  grain  boundaries)  can  be  equivalently  determined, 2*  with  the 
condition  that  in  the  direction  of  the  external  surface  the  surface  energy  vec¬ 
tor  has  magnitude  (ysi  -  ys2  +  Yb/2)»  where  Ysi  =  surface  energy  of  the  parti - 
cle,  ys2  =  surface  energy  of  the  adjoining  particle,  and  Yp  =  common  grain 
boundary  energy. 
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APPENDIX  2 

PARTICLE  ENERGY  FOR  RING  AND  POLYHEDRON  ARRAYS  OF  IDENTICAL  SPHERES 

The  energy  of  particles  in  arrays  formed  by  placing  spheres  in  a  ring 
and  at  the  vertices  of  regular  polyhedra  have  been  calculated  using  the  approach 
outlined  in  Section  2.1,  and  the  four  functions  shown  in  Table  1.  As  outlined 
in  Section  2.1  the  energy  per  particle  is  given  by  E  =  y$[As  +  2A^  cos  (<|*  /2)], 
where  the  surface  area  of  a  particle  (A$),  area  of  the  grain  boundary  (Ab),  and 
particle  radius  (r),  have  been  listed  in  Table  3  and  4  for  the  pertinent  arrays. 

The  function  f-^(z)  is  used  in  defining  the  volume  of  a  spherical  seg¬ 
ment,  illustrated  in  Fig.  Al(a),  with  a  base  radius  of  r  sin  (z/2).  The  func¬ 
tion  f2(x,y,z)  helps  define  the  volume,  while  function  f3(x,y,z)  is  used  to 
define  the  surface  area,  of  the  shaded  segment  shown  in  Fig.  Al(a).^  The  angles 
x,y  and  z,  shown  in  Fig.  Al(a),  will  be  related  to  specific  arrays  in  the 
following  paragraphs. 

Function  f 4 ( 4*,  np,  $,  e)  has  been  used  to  define  the  volume  of  the 
pyramid  with  either  a  pentagon,  square  or  triangular  base  (depending  on  poly¬ 
hedra),  as  shown  in  Fig.  Al(b).  For  polyhedra,  it  is  necessary  to  determine 
this  volume  for  array  configurations  in  which  the  pore  does  not  exist.  The 
angles  9  and  4*  used  in  this  function  are  schematically  shown  in  Fig.  A2  for  a 
cubic  polyhedron;  values  are  listed  by  polyhedra  in  Table  2.  Also  listed  in 
Table  2  is  np,  the  number  of  touching  particles  surrounding  each  particle  in  the 
polyhedron. 

Ring  Array  of  Spheres 

When  a  pore  exists  within  the  ring  array  of  spheres  (4 ,  +  2u/n  <  it), 
the  particle  radius  (r),  surface  area  (As)  and  grain  boundary  area  (Ab)  are 
simply  stated  as  a  function  of  the  contact  angle  (4*),  as  written  in  Table  3. 

For  conditions  where  the  pore  does  not  exist  (4.  +  2u/n  >  u)  the 
geometry  becomes  more  complicated  because  all  particles  penetrate  one  another. 
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The  particle's  radius  (r)  is  calculated  using  functions  f}(z)  and 
f2(x»y>z)»  and  the  particle's  surface  area  (A$)  is  determined  using  the  function 
f3(x,y,z).  For  ring  arrays  the  three  angles  are  given  by: 

x  =  y  =  2  arccos  ( - - - -  ]  , 

ytan  tan  ^  J 


y  =  2n/n  ,  (Al) 

and 

z  = 

The  grain  boundary  area  (Ap)  is  calculated  using  the  contact  angle  (4),  the 
angle  y  (given  by  Eq.  (Al)),  and  the  particle's  radius  (r). 

Polyhedron  Arrays 

For  polyhedra  the  calculations  of  A$,  Ab,  and  r  are  best  divided  into 
different  ranges  of  the  contact  angle  (4>)»  which  describes  different  configura¬ 
tions  of  the  polyhedron.  When  <  u  -  2u/n^  the  pore  within  the  polyhedron  is 
connected  to  the  external  surface  by  a  number  of  channels  equal  to  the  number  of 
faces  on  the  polyhedra.  The  number  of  particles  coordinating  a  channel  n^, 
listed  in  Table  2,  is  equal  to  the  number  of  particles  on  a  polyhedron  face. 
Channels  close  when  4)  =  it  -  2u/n^  creating  an  isolated  pore,  and  the  pore 
disappears  when  4.  =  n  -  2e  (e  are  listed  in  Table  2)  resulting  in  a  fully  dense 
confi gurati on. 

For  the  open  pore  configuration  (4»  <  u  -  2it/n^)  the  particle  radius 
(r),  listed  in  Table  4,  is  calculated  using  f^(z  =  4;),  and  np.  Particle  surface 
area  (A$)  and  grain  boundary  area  (Ab),  listed  in  Table  4,  are  simply  stated  as 
a  function  of  r,  np,  and  4>. 

The  closed  pore  configuration  occurs  for  contact  angles  within  the 
range  u  -  2n/r\f  <  4>  <  u  -  2e.  Particle  radius  (r),  given  in  Table  4,  is 
calculated  using  function  f}  (z  =  4>)  and  f2  (x  =  Y»  y  =  2u/ny,  z  =  4.) ,  where 
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Y  =  2  arccos  ( — s^n  ^ - -)  .  (A2) 

tan  9  tan  ^ 

Particle  surface  area  (Ag)  has  been  calculated  using  function  f$  (x  =  y» 

y  =  2-rt/nf,  z  =  <j;)  and  np.  The  grain  boundary  area  (Ab)  is  a  simple  function  of 

previously  described  variables. 

Geometrical  calculations  become  more  tedious  with  pore  closure 
(<j»  >  it  -  29),  but  generally  the  equations  for  the  particle  radius  (r)  and 
surface  area  (As)  fall  within  three  regimes  dependent  upon  the  contact  angle.  A 
single  equation  describes  the  grain  boundary  area  per  particle  (Ab)  for  the 
three  regimes  (see  Table  A). 

For  the  three  different  regimes  r,  As,  and  Ab  is  calculated  using 
functions  f^(z  =  4>‘ ) ,  fj(z  =  c|>),  f2(x  =  p,  y  =  t),  z  =  4>),  f2(x  =  2u/np, 

y  =  V,  z  =  (];'),  f3(x  =  p,  y  =  t),  z  =  <l») ,  f3(x  =  2u/np ,  y  =  V ,  z  =  <j/ )  and 

f4  (9,  <i>,  np ,  <t>).  The  functional  variables  are  listed  below: 


4»*  =  ±2  arccos 


4>  sin2  C 
7  sin  9 


cos  C 


sin' 


A* 

7 


.  2  2 
sin  <(>  cos 


tan 


(A3) 


(Angles  <j>,  9,  C  are  listed  by  polyhedron  in  Table  2, 
the  condition  cos  (V2)  <  sin  9/tan  C. ) 


2  arctan 


1 


.  <|>  it 

tan  iS—  cos  — 
2  n_ 


t)  =  it  +  V  -  29 


(except  for  the  tetrahedron  when  <j <  <  *r,  then  g  =  n  - 


with  the  terms  negative  for 

(AA) 

(A5) 

V  +  2e) 


p  =  2  arccos 


( tan  -J?  tan 


(A6) 
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Table  4 


Polyhedra 


Open  Pore 

2n 

4,  <  u  -  — 


Closed  Pore 


7t-£i<4,<u-2e 

nf 
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Table  4 
(Continued) 


Polyhedra 


it  .  .  <i>  .  /sin  0-, 

7  -  Q  <f<  arccos 


ri  ^4^9*  4>»  Hp,  <t>)  +  -  itf  ^ ( 4»* )  -  ^2 


(P»  n.  <1») 


R2  ^(1  +  cos  f)  -  ^  (f3(p,  th  4>)  -  f3(^,  V,  4-‘))j 

R2  Tsin2  t  (tc  +  2 4>  -  y  +  sin  y)  +  cos2  J  iilL-liii] 

*  I  £  tan  0  I 


Polyhedra 


!  > arccos 
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Abstract 


It  is  shown  that  grain  boundary  motion  in  a  sintering  array  of  parti¬ 
cles  can  be  limited  by  geometry.  Despite  its  curvature,  the  grain  boundary 
can  encounter  an  energy  barrier  if  it  were  to  move  through  the  smaller  grain. 
This  energy  barrier  decreases  to  2ero  when  the  radius  ratio  (R)  between  adja¬ 
cent  grains  approaches  a  critical  value,  R  -  Rc  =  -(cos  (u»e) )~ 1 ,  where  is 
the  dihedral  angle  relating  surface  and  grain  boundary  energies.  The  change 
in  radius  ratio  (R  -  Rc)  occurs  by  interparticle  mass  transport,  driven  by  the 
differential  surface  curvature  of  the  adjacent  grains.  This  result  suggests 
that  grain  growth  within  a  sintering  array  primarily  occurs  by  interparticle 
mass  transport  making  it  distinctly  different  from  grain  growth  in  the  fully 
dense  material . 

More  significantly,  it  is  shown  that  when  a  grain  disappears  from 
between  other  grains,  the  driving  force  for  mass  transport  to  the  newly  formed 
region  between  the  remaining  grains  is  reinitiated,  i.e.,  the  angle  formed 
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where  the  grain  boundary  meets  the  surfaces  is  smaller  than  the  equilibrium 
(dihedral)  angle.  That  is,  grain  growth  via  interparticle  mass  transport  and 
boundary  motion  can  continually  drive  sintering  and  array  shrinkage.  It  is 
shown  that  the  interparticle  mass  transport  phenomenon  causing  grains  to 
disappear  and  their  neighbors  to  grow  does  not  lead  to  array  shrinkage  when  41 
>  130°,  but  can  lead  to  significant  shrinkage  via  the  sintering  of  remaining 
neighbors.  Thus,  array  shrinkage  and  densif ication  is  predictably  related  to 
grain  growth.  It  is  shown  that  the  relative  density  of  a  ring  of  spheres  that 
undergoes  cyclic  sintering  and  grain  growth  is  nearly  linearly  related  to 
grain  size.  Conceptual  relations  between  density  and  grain  growth  for  real 
powder  compacts  are  discussed  and  related  to  the  function  that  describes  the 
compact's  distribution  of  pore  coordination  numbers.  Finally,  since  interpar¬ 
ticle  mass  transport  producing  grain  growth  can  be  the  rate-limiting  step  for 
densif ication,  the  theory  and  related  experiments  associated  with  densifica- 
tion  kinetics  should  be  reevaluated. 


1.0  INTRODUCTION 
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Grain  growth  in  dense  materials  is  associated  with  grain  boundary 
curvature.  Due  to  curvature,  the  chemical  potential  is  greater  on  one  side, 
causing  atoms  or  ions  to  hop  across.  Because  hopping  (diffusion)  distances  are 
very  small,  boundary  motion  can  be  rapid.  Grain  growth  also  occurs  during  the 
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densif ication  of  a  powder  compact.  Greskovich  and  Lay  have  observed  the  growth 
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of  grains  in  A1 2O3  powder  compacts  with  initial  bulk  densities  as  low  as  30%  of 
theoretical,  e.g.,  upon  densif ication  from  30  to  40%  of  theoretical,  grains  grow 
from  0.3  um  to  0.5  um.  Their  observations  indicated  that  boundaries  will  only 

O 

move  after  necks  have  fully  formed  between  the  sintering  particles.  Gupta's 
analysis  of  literature  data  for  Cu,  A1 2O3 ,  BeO,  and  ZnO  showed  that  grain  size 
increased  very  slowly  during  most  of  the  densification  process,  but  increased 
very  rapidly  once  the  compact's  density  was  >  90  %  of  theoretical.  He  also 
stated  that  within  the  regime  of  slow  grain  growth,  the  grain  size  is  linearly 
proportional  to  density  and  independent  of  the  sintering  temperature.  Bruch 
also  showed  that  the  grain  size  was  only  dependent  on  density  and  independent  on 
temperature.  These  combined  observations  strongly  suggest  that  (a)  grain  growth 
occurs  after  a  substantial  amount  of  mass  transport  has  occurred  to  form  necks, 
(b)  grain  growth  within  the  partially  sintered  body  appears  not  to  be  thermally 
activated,  but  depends  only  on  density,  and  (c)  grain  growth  and  densification 
are,  in  some  way,  intimately  related. 

In  a  previous  paper,  Kellett  and  Lange4  calculated  the  free  energy  and 
shrinkage  of  symmetric  rings  and/or  regular  polyhedra  of  identical  particles 
(cylinders  and  spheres)  that  were  allowed  to  develop  an  equilibrium  config¬ 
uration  by  sintering.  Each  array  contains  a  single  pore.  Configurational 
changes  occurred  when  mass  was  transported  to  the  contact  region  increased  the 
contact  area  (grain  boundary)  between  the  initially  touching  particles.  To 
allow  tractable  analytical  calculations,  it  was  assumed  that  the  curvature  of 
the  particle’s  surface  was  independent  of  surface  position.  That  is,  the  sur¬ 
faces  of  spheres  remained  spherical,  but  increased  their  radius  to  conserve  mass 
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as  the  particle  penetrated  one  another.  It  was  also  assumed  that  each  particle 
would  undergo  identical  shape  changes.  This  assumption  implicitly  requires  no 
net  transport  of  mass  between  particles  and  that  grain  boundaries  remain  planar. 
The  free  energy  per  particle  was  determined  as  a  function  of  array's  configura¬ 
tional  change  by  summing  the  energies  associated  with  its  surface  and  grain 
boundaries.  Young's  equation  was  used  to  equate  the  grain  boundary  energy  per 
unit  area  to  the  particle's  external  surface  energy  per  unit  area.  The  angle, 

4),  defined  by  the  tangents  to  the  free  surface  where  the  grain  boundary  meets 
the  free  surface,  was  used  to  define  the  configurational  change  (ip  =  0  when  the 
particles  are  just  touching  and  increases  as  the  spheres  penetrate  one  another). 

The  results  pertinent  to  this  paper  are:  the  free  energy  of  the  arrays 
depends  on  conf iguration  (defined  by  ip)  and  exhibits  a  minimum  when  the  contact 
angle  equals  the  dihedral  angle  (ip  =  ipe) ;  the  dihedral  angle  is  defined  by  the 
surface  to  grain  boundary  energy  ratio  through  Young’s  equation.  When  the 
array's  coordination  number  exceeds  a  critical  value  (n  >  nc)  the  equilibrium 
configurations  still  contain  a  pore,  i.e.  the  pore  shrinks  to  an  equilibrium 
size.  When  n  <  nc,  the  pore  disappears  as  the  array  shrinks  to  its  equilibrium 
configuration;  nc  increases  with  the  dihedral  angle.  The  same  result  concerning 
pore  stability/disappearance  is  reached  when  an  isolated  pore  in  a  finite  poly¬ 
crystalline  body  is  analyzed  by  a  similar  method,  e.g.,  the  isolated  pore  will 
be  stable  when  the  surface  curvature  of  the  grains  surrounding  the  pore  is 
identical  to  the  surface  curvature  of  grains  on  the  external  surface.5 

5 

In  another  paper,  it  was  hypothesized  that  grain  growth  would  decrease 
the  coordination  number  of  stable  pores  to  the  condition  where  they  would  be 
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thermodynamically  unstable  and  disappear  (i.e.,  where  n  <  nc).  This  hypothesis 
was  tested  by  studying  the  grain  growth  required  to  eliminate  mono-sized  spheri¬ 
cal  voids  in  an  otherwise  dense  body.  The  spherical  voids  were  produced  by 
incorporating  plastic  spheres  into  the  powder  compact  prior  to  sintering.  It 
was  observed  that  the  number  density  of  spherical  voids  decreased  with  grain 
growth,  whereas  the  voids  that  remained  after  a  given  period  of  grain  growth  did 
not  change  in  size.  These  observations  are  consistent  with  the  initial  hypoth¬ 
esis,  suggesting  that  grain  growth  was  helpful  in  eliminating  thermodynamically 
stable  voids. 

The  object  of  this  paper  is  to  examine  the  conditions  for  grain  growth 
in  a  powder  compact  and  to  relate  this  grain  growth  phenomena  to  the  compact's 
sintering  behavior.  Arrays  of  spheres  will  be  used  to  examine  the  conditions 
for  grain  growth,  and  densif ication  produced  by  grain  growth  in  a  similar  manner 
used  previously.5  It  will  be  shown  that  (I)  grain  growth  in  a  powder  compact  is 
not  a  simple  matter  of  boundary  motion,  (2)  grain  growth  reinitiates  sintering 
and  therefore  densif ication,  and  (3)  the  amount  of  grain  growth  required  to 
produce  a  dense  body  is  related  to  the  pore  coordination  number  distribution 
within  the  initial  powder  compact. 

2.0  CONFIGURATIONAL  CHANGES  FOR  GRAIN  BOUNDARY  MOTION 

In  the  current  work,  the  initial  mass  of  adjoining  particles  will  be 
different,  producing  differential  free  energies  and  the  development  of  nonplanar 
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grain  boundaries.  The  object  of  this  section  is  to  examine  the  configurational 
change  required  to  produce  an  incremental  decrease  in  free  energy  for  an  incre¬ 
mental  displacement  of  the  grain  boundary,  i.e.,  the  configuration  where  grain 
boundary  motion  will  continually  decrease  the  free  energy  of  the  system. 

In  temperature  regimes  for  mass  transport,  the  net  negative  curvature 
that  develops  at  the  contact  position  between  touching  particles  results  in  a 
large  chemical  potential  difference,  relative  to  other  surface  regions,  and  thus 
a  large  driving  force  for  material  transport  to  the  region  where  the  particles 
touch  and  a  grain  boundary  develops.  Mass  transported  to  the  contact  (or  neck) 
region  can  primarily  come  from  one  of  the  adjoining  particles  if  their  initial 
free  energies  (i.e.,  surface  curvatures)  are  very  different,  or  from  both  par¬ 
ticles.  When  the  size  of  the  adjoining  particles  are  different,  their  differ¬ 
ential  surface  curvature  can  also  result  in  interparticle  mass  transport.  But, 
because  the  driving  force  causing  mass  transport  to  the  neck  region  is  much 
greater,  it  might  be  assumed  little  net  transport  between  the  particles  will 
occur  until  the  driving  force  to  the  neck  region  is  dissipated  through  configu¬ 
rational  change.  For  this  reason,  although  intraparticle  transport  producing 
neck  growth  (sintering)  and  interparticle  transport  producing  grain  growth  may 
be  concurrent  within  a  powder  compact,  one  might  expect  that  when  the  driving 
force  exists,  mass  transport  to  the  neck  region  will  dominate  relative  to  mass 
transport  between  particles. 

For  all  arrays  considered,  it  is  assumed  that  the  curvature  of  each 
particle's  surface  is  independent  of  surface  position,  e.g.,  spherical  surfaces 
remain  spherical.  With  this  assumption,  touching  spheres  within  the  array 
change  their  configuration  by  penetrating  one  another.  As  they  penetrate,  their 
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radii  change  either  to  conserve  their  mass  (condition  of  intraparticle  mass 
transport)  or  change  their  mass  as  one  particle  grows  at  the  expense  of  a  neigh¬ 
boring  particle  (interparticle  mass  transport).  Although  this  type  of  configu¬ 
rational  change  is  not  theoretically  expected  or  experimentally  observed  during 
the  initial  stages  of  neck  growth  (changes  are  much  more  rapid  closer  to  the 
neck  region),  it  does  approximate  the  real  situation  and  allows  tractable 
analytical  calculations. 

2.1  Two  Particles 

Figure  1  illustrates  the  expected  configurational  change  for  two 
touching  spheres  with  the  assumptions  that  mass  is  only  transported  to  the 
neck  region.  As  detailed  elsewhere,4  the  contact  angle,  n,  is  a  convenient 
parameter  to  describe  the  configurational  change.  When  the  spheres  are  just 
touching,  it>  =  0.  A  previous  analysis  has  shown  that  the  particles  will  inter¬ 
penetrate  until  the  contact  angle  equals  the  dihedral  angle,  i.e.,  i>  =  *  . 
Because  each  particle  conserves  its  mass,  each  increases  its  radius  of  curva¬ 
ture.  Using  the  geometric  parameters  shown  in  Fig.  2a,  the  radius  of  each 
particle  (r^)*  norma1ized  by  its  initial  radius  (r^  g)  can  be  determined 
using  the  following  equations: 


R1  =  ri/r i 1  =  [1-0.25(1  -  cose1)2(2  +  cose^) 

+  2(^t)3cos3(4./2)(1  -  cos»)2(2  +  cos*)]"173  (1) 
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R2  =  r2/r i 2  =  l1  ”  °*25  '  cose2)2(2  +  cose2) 


-  — 7— ,  cos3(*/2)(1  -  cos*)2(2  +  cos,)]’173  (2) 

(R-ir 


where  ej  =  arctan  (Rcj^Ti> 

@2  =  i*  ~  0  ^ 

<t>  =  e i  -  i \>/2. 

R  defines  the  radius  ratio  of  the  two  particles  (R  =  r 2/ri)  during  their 
configurational  change.  The  initial  radius  ratio  (R^  =  2/r i 1 )  can  be 

calculated  using 

Ri  =  R  (R1/R2).  (3) 

The  radius  of  curvature  for  the  grain  boundary  is  given  by 

rb  =  2r1r2/(r2  -  rj)  cos(u>/2)  =  [  2r2/ (R  -  1)]  cos  («i»/2 ) .  (4) 

Before  we  examine  the  configurational  changes  produced  by  mass  transport 
between  the  particles,  let  us  examine  the  conditions  for  grain  boundary  motion 
that  can  produce  a  continuous  decrease  in  the  system's  free  energy. 
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Figure  1  shows  that  one  of  two  configurations  can  be  reached  when 
ip  =  ipe.  The  upper  configuration  (c)  occurs  when  R  <  -(cos(ipe))-^.  For  this 
configuration,  the  two  particles  could  lower  their  free  energy  by  allowing  the 
grain  boundary  to  become  planar,  but  if  the  boundary  were  to  move  through  the 
smaller  grain,  it  would  have  to  increase  its  area  and  thus  the  free  energy  of 
the  system.  That  is,  an  energy  barrier  exists  for  grain  boundary  motion  when 
particles  develop  configuration  (c)  during  sintering.  Configuration  (d), 
occurs  when  R  >  -(cos(ipe))"* .  Here  it  is  obvious  that  the  grain  boundary  can 
move  through  the  smaller  particle  and  decrease  its  surface  area  (and  the  free 
energy  of  the  system)  either  with  or  without  changing  its  curvature.  From 
this  configuration,  grain  boundary  motion  will  cause  a  continuous  decrease  in 
free  energy.  This  configuration  occurs  when  the  size  ratio  of  the  particles, 
developed  during  the  configurational  change,  satisfies  the  relation 

R  =  Rc  =  -  (cos  (ip) ) -1 ,  (5) 

where  Rc  is  the  critical  radius  ratio  which  allows  the  grain  boundary  to  move 
through  the  smaller  grain  and  continually  decrease  its  surface  area. 

Figure  2b  details  the  geometry  for  the  condition  which  just  allows 
the  grain  boundary  to  move  through  the  smaller  grain  without  encountering  an 
energy  barrier.  Note  that  for  this  special  configuration,  the  center  of 
curvature  of  the  smaller  grain  lies  on  the  imaginary  plane  (broken  line  in 
Fig.  2b)  separating  the  two  grains.  Also,  for  this  configuration  the  radius 
of  curvature  of  the  grain  boundary  is  given  by 
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rg'  =  -r£  (cos»p/cos(4>/2)  ) .  (6) 

Thus,  it  can  concluded  that  when  the  size  ratio  of  the  two  particles  becomes 
greater  than  the  critical  ratio,  i.e.,  R  >  Rc  =  -(cos(4ie))“* ,  the  grain 
boundary  can  move  through  the  small  grain  and  continuously  decrease  the  free 
energy  of  the  system  without  mass  transport  between  the  two  particles.  Fig¬ 
ure  3  illustrates  the  critical  size  ratio  (expressed  as  the  initial  size 
ratio,  R.j)  as  a  function  of  the  dihedral  angle  (4^).  For  dihedral  angles 
ranging  from  100-150°,  the  initial  particle  size  ratio,  which  can  cause  the 
grain  boundary  to  spontaneously  move  through  the  small  grain,  decreases  from 
6.6  to  1.1. 

Figure  1(c)  shows  that  the  differential  free  energy  produced  by  the 
different  particle  curvatures  can  cause  mass  to  be  transported  from  the 
smaller  to  the  larger  particle.  During  interparticle  mass  transport,  the 
radius  of  the  larger  particle  increases  as  that  of  the  smaller  decreases. 
Since  material  has  already  been  transported  to  the  neck  region,  the  dihedral 
angle  remains  constant.  The  increase  in  the  radius  ratio  will  cause  R  -►  Rc 
shown  by  step  (d).  Fig.  1,  and  the  condition  when  the  grain  boundary  can  move 
through  the  smaller  grain  while  decreasing  the  free  energy  of  the  system. 

Once  the  grain  boundary  has  disappeared  leaving  the  unusually  shaped  single 
grain  (step  (e).  Fig.  1)  it  can  further  decrease  its  free  energy  by 
spheridizing. 
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2.2  Three  Particles 

Three  touching  spheres  can  be  arranged  such  that  their  centers  are 
either  col  inear  or  noncoliner.  Configurational  changes  will  be  considered  for 
both  cases  since  they  produce  simply  visualized  results  pertinant  to  arrays 
that  enclose  void  space. 

2.2.1  Three  Col  inear  Particles 

Two  of  the  three  particles  are  assumed  to  be  identical.  The  case 
where  the  third  particle  is  larger  and  sandwiched  between  the  two  smaller 
particles  is  nearly  identical  to  the  case  for  the  two  particles  considered 
above.  The  radius  of  the  sandwiched  particle  increases  more  rapidly  with 
changing  because  of  the  two  interpenetrating  smaller  particles.  The  radius 
ratio  where  the  boundaries  can  move  through  the  small  grains  will  still  be 
defined  by  Eq.  (5),  but  this  will  be  developed  by  a  somewhat  smaller  initial 
radius  ratio  than  for  the  case  of  the  two  spheres. 

The  most  interesting  case  is  where  a  smaller  particle  is  sandwiched 
by  two  larger  particles  schematically  shown  in  Fig.  4.  We  will  only  detail 
the  case  where  R  <  Rc  =  -(cos  (4<e))"^.  Configurational  changes  from  where  the 
particles  are  touching  (hj  =  0)  to  where  i|>  -►  ipe  are  shown  in  steps  (a)  to  (b) 
in  Fig.  4.  This  configurational  change,  developed  by  mass  transport  to  the 
neck  region  is  similar  to  that  for  the  two  particles  except  that  since  the 
smaller  particle  is  penetrated  by  two  particles  its  radius  ratio  increases 
more  rapidly  with  i>  (while  conserving  its  mass).  The  radius  of  the  grain 
boundaries  is  given  by  Eq.  (4).  Step  (b)  represents  the  configuration  where 
the  driving  force  for  mass  transport  to  the  neck  region  has  diminished  to 
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zero,  i.e.,  the  equilibrium  contact  angle  (iie)  has  been  achieved.  In  this 
configuration,  an  energy  barrier  exists  for  grain  boundary  motion. 

Interparticle  diffusion,  driven  by  the  differential  surface  curvature 
between  the  smaller  and  larger  grains,  will  increase  the  radius  ratio  and 
result  in  the  configuration  shown  by  step  (c).  Fig.  4.  This  configuration 
develops  when  the  radius  ratio  increases  to  R  =  Rc  =  - (cos (^e) ^ .  i.e.,  the 
condition  found  for  two  particles  where  the  grain  boundary  could  move  without 
encountering  an  energy  barrier.  It  should  be  noted  that  when  R  =  Rc,  the  two 
larger  particles  just  touch  one  another  but  are  still  separated  by  the  smaller 
grain  (shown  by  broken  lines).  It  can  be  shown  (see  Fig.  6,  B  =  n)  that  the 
contact  angle  between  the  two  larger  particles  is  given  by 
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Both  grain  boundaries  can  move  without  encountering  an  energy  barrier  to 
produce  a  single  boundary  that  separates  the  two  larger  grains,  as  shown  by 
step  (c).  Fig.  4.  Since  the  contact  angle  developed  between  the  larger  grains 
by  interparticle  mass  transport  is  less  than  the  dihedral  angle,  viz., 
i|i  =  2ii>e  -u  <  i|ie,  a  new  driving  force  exists  for  mass  transport  to  the  neck 
region  between  the  larger  particles.  Further  mass  transport  to  the  neck 
region  will  cause  4>  -►  i(ie  and  the  development  of  the  stable  configuration  shown 
by  step  (d),  Fig.  4. 

This  simplistic  analysis  illustrates  that  the  disappearance  of  the 
smaller  grain  from  between  the  larger  grains  produces  a  condition  for  mass 
transport  to  the  neck  region  between  the  larger  grains.  That  is,  grain  growth 
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produced  by  interparticle  mass  transport  can  reestablish  the  stronger  driving 
force  for  sintering. 

2.2.2  Three  Noncolinear  Particles 

Figure  5  illustrates  the  configurational  change  for  the  case  where 
the  three  particles  are  noncolinear  and  the  smaller  particle  separates  the 
two,  larger  identical  particles;  e  defines  the  angle  between  the  centers. 
Configuration  (b)  is  produced  as  the  contact  angle  increases  and  the  distance 
between  the  larger  and  smaller  particles  decreases  due  to  mass  transport  to 
the  neck  regions.  As  the  contact  angle  continues  to  increase,  a  special 
configuration  can  be  produced  where  the  two  larger  spheres  just  touch  one 
another  at  one  position  (configuration  c).  It  can  be  shown  that  this  config¬ 
uration  will  develop  whenever  8  <  2  arctan(R  simi*e/(R  cosit>e  +  1)),  or  for  the 
case  where  R  =  1,  when  8  <  ipe .  If  this  configuration  is  not  produced  by  mass 
transport  to  the  neck  region  (i.e;  6  >  i^e)  it  will  be  produced  eventually  by 
interparticle  mass  transport.  When  the  particles  just  touch,  it  can  be  shown 
(see  geometric  of  Fig.  6)  that  the  contact  angle  formed  between  them  can  be 
expressed  as 

si  rbf 

=2it»-it  =  8-  n  +  2  arcsin  ( — (8) 
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Since  n»‘  <  i,  a  greater  driving  force  for  mass  transport  will  exist  at  the 
neck  region  where  the  two  larger  particles  touch  relative  to  where  the  larger 
particles  are  in  contact  with  the  smaller  particle.  It  should  be  pointed  out 
that  the  same  situation  is  developed  when  all  particles  are  identical,  R ^  =  1 
and  e  <  4>e  - 

Interpenetration  of  the  two  larger  particles,  which  initiates  at 
configuration  (c)  and  ends  with  configuration  (e),  causes  e  to  decrease,  i.e., 
the  two  larger  particles  rotate  toward  each  other  about  the  smaller  parti¬ 
cle.*  (In  addition,  if  the  particles  are  cylinders,  a  portion  of  the  bounda¬ 
ries  between  the  larger  and  smaller  cylinders  can  move  to  form  a  planar  bound¬ 
ary  between  a  portion  of  the  larger  cylinders  without  encountering  an  energy 
barrier;  a  configuration  similar  to  that  shown  by  (d)  is  formed.  This  cannot 
occur  for  spherical  particles,  and  thus  the  details  will  not  be  presented.) 

Configuration  (e)  develops  when  -  vpe  at  all  contact  positions.  If 
the  condition  expressed  by  Eq.  (5)  was  not  met  before  tp  -  ipe,  then  inter¬ 
particle  mass  transport  will  take  place,  driven  by  the  different  particle 
curvatures,  to  produce  configuration  (e)  where  the  two  curved  boundaries  can 


*  As  shown  previously,4  an  infinite  linear  array  of  identical  particles  will 
form  an  equilibrium  configuration  where  the  surface  curvature  of  every 
grain  is  identical  to  one  another.  If,  on  the  other  hand,  the  same  linear 
array  is  kinked  such  that  the  angle  s  between  the  three  particles  that  form 
the  kink  satifies  the  relation  b  <  *-,  then  it  can  be  shown  that  the  three 
particles  will  develop  surface  curvature  different  than  the  rest.  This 
condition  arises  because  the  three  penetrate  one  another,  whereas,  all 
others  only  penetrate  their  neighbors.  Thus,  even  when  the  particles  are 
initially  identical  and  conserve  their  mass  during  neck  growth,  differen¬ 
tial  surface  curvature  can  develop  leading  to  interparticle  mass  transport 
and  grain  growth. 
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move  together  without  encountering  an  energy  barrier.  Final  mass  transport 
can  result  in  the  same  configuration  shown  as  (d)  in  Fig.  4,  i.e.,  two 
identical,  penetrating  spheres  separated  by  a  planar  boundary. 


V. 

/  V  /  y 
iWj 


3.0  GRAIN  GROWTH  AND  ARRAY  SHRINKAGE 


<•.  •  -  •'1 


The  above  analysis  showed  that  it  is  energetically  favorable  for 
grain  boundaries  to  move  across  sintering  particles  only  (1)  after  substantial 
neck  growth  and  (2)  when  the  radius  ratio  of  adjacent  grains  exceeds  a  criti¬ 
cal  value  (see  Eq.  (5)).  The  first  of  these  two  conditions  was  implied  by 
Greskovich  and  Lay.1  It  was  also  shown  that  when  the  initial  particles  are 
somewhat  similar  in  size  (see  Fig.  3),  interparticle  mass  transport  is  first 
required  to  develop  the  critical  radius  ratio  for  boundary  motion.  More 
significantly,  it  was  shown  that  when  smaller  grains  disappear  between  larger 
grains,  a  new  (or  increased)  driving  force  is  established  for  neck  growth 
between  the  larger  grains.  That  is,  the  disappearance  of  smaller  grains 
produced  by  a  combination  of  interparticle  mass  transport  and  boundary  motion 
can  either  reestablish  or  enhance  the  driving  force  for  mass  transport  to  the 
neck  region,  viz.,  grain  growth  can  enhance  the  sinterabil ity  of  a  powder 
compact.  The  object  of  this  section  will  be  to  establish  a  relationship 
between  grain  growth  and  shrinkage  of  simple  particle  arrays. 
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3.1  Linear  Arrays 

Assume  that  the  configurations  in  Fig.  4  are  repeated  to  form  a 
linear  array  with  two  separating  sets  of  spheres,  with  initial  radii  r^  and 
r^*  Configuration  (b),  developed  by  mass  transport  to  the  neck  regions,  is 
assumed  to  produce  a  dihedral  angle  where  the  grain  boundaries  intersect  the 
surface,  i.e.,  R  <  Rc.  This  requires  that  the  initial  radius  ratio 
(Ri  =  r^/r^)  is  less  than  a  given  size  (similar  to  that  illustrated  in  Fig.  3; 
for  the  linear  array,  each  particle  penetrates  two  others  which  changes  the 
relation  between  the  initial  and  final  radius  under  conditions  of  mass  conser¬ 
vation  from  that  for  the  two  penetrating  spheres  detailed  above).  Configura¬ 
tion  (c)  develops  by  interparticle  mass  transport  causing  the  larger  particles 
(defined  by  subscript  2)  to  increase  their  mass  at  the  expense  of  diminishing 
the  mass  of  the  smaller  particles.  Interparticle  mass  transport  causes  R  -  Rc 
where  the  two  boundaries  move  together  to  form  a  single  boundary  between  the 
two  touching,  larger  particles  as  discussed  above  for  the  case  of  three 
colinear  particles.  The  contact  angle,  =  2ipe  -  it,  developed  between  the  two 
larger  particles  in  configuration  (c)  is  given  by  Eq.  (7).  Since  <  *e,  mass 
transport  to  the  neck  regions  is  reinitiated  to  develop  the  dihedral  angle 
between  the  larger  particles  as  shown  in  configuration  (d). 

The  shrinkage  of  the  linear  array  can  be  determined  for  each  of  the 
configurational  steps.  Since  engineering  strains  are  used  to  characterize  the 
shrinkage  of  powder  compacts  during  densif ication,  the  shrinkage  strain  of  the 
linear  array  will  be  determined  by 

e  =  *  ^b,c,d^a* 
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where  the  initial  center  to  center  distance  between  the  touching  larger 
spheres  (configuration  a)  is  given  by  la  =  r^  +  ri2  =  ru(Ri  +  1)  and  lb  , 
lc,  and  ld  are  the  distances  between  the  larger  grains  for  configurations  (b) 
through  (d)  and  are  given  by  the  following  equations: 


lb  =  r^R^lcose^  +  Rcos8  2  ] 

(10a) 

!c  =  rn(l  +  R3)1/3  F(2^e  -  *)  s i nii>e 

(10b) 

ld  -  rn(  1  +  R3) 1/3  F(*e)  cos Ue/2) 

(10c) 

where 

F(z)  =  [0.5  cos(§)(3  -  cos2(f)))-1/3  , 


R1  =  [1  -  0.5(1  -  C0S81)2(2  +  COSOj) 

+  4(|— j)3  cos3(^|)(1  -  cos*)2(2  +  cos,)]-1/3 


R2  =  [1  -0.5(1  -  cose2)2(2  +  cose2) 

- — ,  cos3(^)(l  -  cos<t>) 2 (2  +  cos  o)]*1/3 

(R  -  l)3  2 

R  si  n«p 

ei  =  arctan?rros»e~ri1 
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02  =  V  61 


6^  -  t-g/2  , 


and  R  =  R^ ( R2/R  x) 


Figure  7  illustrates  the  shrinkage  strains  (R^  =  1)  for  the  three  configura¬ 
tions,  related  to  the  initial  configuration,  as  a  function  of  the  dihedral 
angle,  i|>e.  As  shown,  the  major  contributions  to  shrinkage  occur  when  mass  is 
transported  to  the  neck  regions,  i.e.,  between  configurations  (a)  and  (b),  and 
(c)  and  (d).  Interparticle  mass  transport,  which  causes  the  smaller  particles 
to  disappear,  also  causes  the  larger  particles  to  grow  and  thus  does  not  pro¬ 
duce  significant  shrinkage  of  the  linear  array.  Although  interparticle  mass 
transport  does  not  directly  contribute  to  shrinkage,  it  does  lead  to  the 
reinitiation  of  neck  growth  which  produces  significant  shrinkage  once  the 
smaller  particles  disappear.  Thus,  it  might  be  concluded  that  grain  growth 
within  the  linear  array,  which  occurs  by  a  combination  of  interparticle  mass 
transport  and  boundary  motion,  reinitiates  the  driving  force  for  sintering  and 
thus,  shrinkage. 

One  can  show  that  if  the  initial  radius  ratio  Is  large  enough  such 
that  R  >  Rc,  configuration  (c)  is  reached  directly  from  (a)  without  inter¬ 
particle  mass  transport.  The  same  conclusions  concerning  shrinkage  are 
obtained  without  interparticle  mass  transport.  It  can  also  be  shown  that  if 
R.|  >  1,  but  R  <  Rc,  the  same  conclusions  are  obtained  as  shown  in  Fig.  8 
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except  that  the  shrinkage  strain  for  each  configurational  change  decreases 
with  increasing  R^. 

The  grain  growth  phenomena  and  its  associated  shrinkage  between 
configurations  (b)  and  (d)  can  be  repetitive.  To  illustrate  this,  let  us 
consider  a  finite  linear  array  of  nQ  nearly  identical  spheres  of  radius  ri , 
i.e.,  Ri  »  1.  The  initial  length  of  the  array  is  L0  =  2rinQ.  Mass  transport 
to  the  contact  regions  will  result  in  conf iguration  (b)  and  a  shrinkage  of  the 
array  to  a  new  length  defined  as  =  2r^  F(^e)  cos(nie/2)  n0,  where  the  function 

F(\^e)  defines  the  new  radius  of  the  penetrating  spheres  in  terms  of  the  initial 
radius,  i.e.,  r  =  r.  F  (^e)  =  r.  [cos(*e/2)(3  -  cos2  (i^e/2) )  ] Now  assume  that 
interparticle  mass  transport  occurs  such  that  every  other  grain  disappears  and 
the  array  develops  configuration  (d).  During  this  process,  the  number  of  grains 
is  reduced  by  a  factor  of  2  or,  conversely,  the  new  grains  increase  their  radius 
by  a  factor  of  2^  (configuration  (c)).  Upon  further  mass  transport  to  the 
contact  regions  (configuration  (c)  to  (d)),  the  array  shrinks  to  a  new  length 
given  by  L2  =  2ri(21^3)  F(^e )  cos(i^e/2)  (nQ/2).  Configuration  (d)  now  reiniti¬ 
ates  a  new  cycle  of  grain  growth  via  interparticle  mass  transport  and  boundary 
motion  as  configuration  (b).  During  each  cycle,  the  number  of  grains  is  reduced 
by  a  factor  of  2,  and  the  array  shrinks  by  the  difference  between  the  mass 
centers  of  configurations  (b)  and  (d).  At  the  end  of  each  cycle  (defined  by  m) 
it  can  be  shown  that  the  grain  size  is  given  by 

Dm  -  2ri(2m/3),  (11) 
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and  the  array  length  by 


W  ■  2""(V2rl)  F(*e>  cos(*,e/2)L0. 


Equation  (12)  can  be  expressed  as  the  shrinkage  strain  of  the  array,  relating 
shrinkage  strain  after  m  cycles  of  grain  growth  to  grain  size: 


‘m+1  =  WLo  -  1  =  2_m(V2ri)  cos(^e/2)  -  1. 


3.2  Closed  Arrays 

Unlike  linear  arrays,  closed  arrays,  viz,  rings  of  cylinder  or 
spheres  and  polyhedra,  contain  a  void  space  that  develops  into  a  pore  upon 
sintering.  The  number  of  particles  coordinating  the  void  (or  pore)  is  defined 
as  the  pore  coordination  number,  n.  As  discussed  above,  it  was  shown4  that 
during  sintering,  the  pore  formed  with  identical  particles  can  either  disap¬ 
pear  (n  s  nc)  or  shrink  to  an  equilibrium  size  (n  >  nc),  where  nc  is  a  criti¬ 
cal  coordination  number  which  depends  on  the  dihedral  angle  (nc  decreases  as 
4>e  increases).  It  has  been  suggested5,6  that  grain  growth  will  lower  the 
coordination  number  of  stable  pores  (n  >  nc),  such  that  it  will  disappear 
(when  n  <  nc).  Using  the  iterative  grain  growth  phenomena  described  above, 
the  relation  between  grain  growth,  shrinkage  (or  densification) ,  and  pore 
disappearance  will  be  detailed  below  for  the  ring  of  spheres. 

Consider  a  symmetric  ring  containing  two  different  sets  of  spheres 
with  nearly  the  same  radius,  r^  (Ri  -  1);  the  total  number  of  spheres  is  nQ. 
The  total  mass  of  material  in  the  ring  is  nn(4/3)  n  r^3  pt,  where  pt  is  the 
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material  density.  The  ring's  volume  can  be  described  by  a  disc  with  a  radius 
that  circumscribes  the  ring  and  a  height  equal  to  the  sphere  diameter;  its 
initial  volume  is  given  by  2n  r^( 1  +  (sin(eQ/2))  )  ,  where  e0  is  the  angle 

which  initially  subtends  one  particle,  i.e.,  e0  =  2n/n0  .  Upon  sintering,  the 
radius  of  the  circle  that  circumscribes  the  array  decreases  to  r^  F ( ipe /2 )  (1  + 
cos(u>e/2)(sin(e0/2))-1)  and  its  volume  becomes  2n  r^  F ( »Pe/2 ) 3  (1  +  cos(4>e/2) 

1  p 

(sin(e0/2))~  )  .  Since  nQ  >  nc,  the  pore  within  the  array  does  not  disappear 
during  this  initial  step  of  sintering. 

Interparticle  mass  transport  can  now  dominate  which  is  assumed  to 
cause  every  other  grain  to  disappear,  as  described  for  the  case  of  the  nonco- 
linear  particles  above.  (If  the  initial  ring  was  not  symmetrical,  the  analy¬ 
sis  for  the  case  of  the  noncol inear  particles  shows  that  torques  would  develop 
during  grain  growth  causing  the  ring  to  distort  from  its  initial  configura¬ 
tion.)  During  the  first  cycle  of  grain  growth,  the  pore's  coordination  number 
is  reduced  by  a  factor  of  two  (n^  =  nQ/2),  the  angle  subtending  each  grain  is 
increased  by  the  same  factor  (i.e.,  e  =  2e0  =  4Ti/n0),  and  the  radius  of  the 
circumscribed  circle  decreases  by  a  factor  of  2^3  (i.e.,  =  21/3  r^  F(4>e/2) 

(1  +  cos(4'e/2) (sin(e^/2))-^).  As  for  the  arrays  discussed  above,  it  is  not 
the  disappearance  of  the  grains,  but  the  reinitiation  of  sintering  that  causes 
the  array  to  shrink  upon  grain  growth.  If  n^  >  nc,  the  pore  will  still  remain 
after  the  first  cycle  of  grain  growth;  subsequent  cycles  of  grain  growth  will 
cause  n  <  nc  and  pore  disappearance. 

Subsequent  cycles  (m)  of  grain  growth  will  cause  the  volume  of  the 
ring  array  to  decrease  in  a  manner  similar  to  that  described  for  the  first 
cycle.  Since  its  mass  remains  constant,  its  density  increases  with  grain 
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growth.  It  can  be  shown  that  the  density  of  the  ring  array  can  be  expressed 
as 

pm+1  =  (2/3) (Dm/2ri)-3  F(*e/2)-3  (1  +  cos(*e/2) (sin(em/2))-1)-2n0pt, 

(14) 

where  Dm/2ri  is  the  grain  size  normalized  to  the  initial  sphere  diameter  (see 
Eq.  (11)),  em  =  (2m)  e0  =  2*/nm,  and  nm  =  (2_m)n0.  The  pore  will  just  close 
when  nm  =  nc  (i.e,  when  em  =  -n  -  4<e) .  The  density  of  the  array  when  the  pore 
just  closes  (chosen  as  the  theoretical  density  of  the  array*)  can  be  expressed 
as 

Pt (array)  =  (2/3)  F(^/2)"3(l  +  cos(*e/2)(sin((i.  -  o^)^))'1)"^^. 

(15) 


Figure  9  illustrates  the  relative  density  of  the  ring  array 
(prtH-l/pt (array))  as  a  grain  size  normalized  by  the  diameter  of  the 

initial  spherical  particles.  As  illustrated,  the  near  linear  relationship 
between  relative  density  and  grain  size  over  much  of  the  densification  range 
is  remarkably  similar  to  the  experimental  relations  reported  by  Gupta.  Also, 


*  Further  grain  growth  will  result  in  the  interpenetration  of  all  grains  and 
the  production  of  one  spherical  grain.  The  one  spherical  grain  only  occupies 
2/3  of  the  disc  volume  used  in  the  array  calculations  and  would  appear  to 
cause  the  density  of  the  array  to  decrease  when  rL,  <  nc.  This  artifact  of  the 
calculation  can  be  avoided  by  either  redefining  tne  array  volume  by  the  sphere 
that  circumscribes  the  array  or  by  choosing  the  array's  theoretical  density  as 
that  expressed  by  Eq.  (15). 
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the  amount  of  grain  growth  that  takes  place  during  densification  is  also 
remarkably  similar  to  the  experimental  data.  Figure  9  also  shows  that  more 
grain  growth  is  required  for  densification  with  decreasing  dihedral  angle. 


DISCUSSION 

It  has  been  shown  that  a  special  geometrical  relation  must  exist 
(R  >  Rc)  between  two  adjoining  grains  before  their  grain  boundary  can  move, 
without  encountering  an  energy  barrier.  This  result  strongly  suggests  that 
the  motion  of  grain  boundaries  will  be  governed  by  geometrical  considerations 
within  a  densifying  powder  compact,  i.e.,  its  relative  density.  Because  a 
given  relative  density  can  be  achieved  by  different  temperature/time 
schedules,  grain  growth  will  appear  to  be  unrelated  to  these  different 
schedules. 

As  shown,  grain  growth  in  a  densifying  powder  compact  is  primarily  a 
result  of  interparticle  mass  transport.  Since  the  diffusion  distance  for 
interparticle  mass  transport  is  some  fraction  of  the  grain  size,  whereas  that 
for  boundary  motion  is  related  to  atomic  distances,  one  would  expect  grain 
growth  kinetics  in  densifying  bodies  to  be  much  slower  than  when  the  same  body 
is  ful ly  dense. 

The  most  significant  observation  of  the  current  work  is  that  when 
grains  disappear  from  between  adjoining  grains  by  a  combination  of  interparti¬ 
cle  mass  transport  and  boundary  motion,  the  driving  force  for  mass  transport 
to  the  neck  regions  is  reinitiated.  That  is,  once  the  initial  particles 
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sinter  together,  further  densification  is  driven  by  grain  growth.  This  obser¬ 
vation  leads  to  the  near  linear  relation  between  density  and  grain  size  for 
the  ring  array  of  particles.  The  following  will  discuss  other  implications 
concerning  densification  and  grain  growth  in  powder  compacts  and  densification 
kinetics  can  be  implied  from  this  observation. 

In  a  previous  paper, 6  it  was  suggested  that  the  thermodynamics  of 
densification  are  controlled  by  the  function  that  describes  the  distribution 
of  pore  coordination  numbers  within  a  powder  compact.  This  function,  V(n), 
describes  the  volume  fraction  of  pores  associated  with  each  and  every  coordi¬ 
nation  number.  The  total  volume  fraction  of  void  phase  within  the  compact  is 
the  sum  of  the  coordination  number  distribution  function  over  all  coordination 
numbers.  For  example,  if  identical  particles  are  packed  with  a  simple  cubic 
arrangement,  all  pores  are  coordinated  by  8  particles.  The  function  V(n)  is 
only  nonzero  when  n  =  8,  viz.,  V(8)  =  1  -  n/6  =  0.48  (the  volume  fraction  of 
voids  within  this  arrangement).  If  this  particle  arrangement  was  sintered  for 
the  case  where  u>e  ^  109.5°  (or  nc  <  8),  all  pores  would  be  thermodynamically 
unstable  and  disappear,  kinetics  permitting.  For  this  case,  pores  do  not  need 
to  be  recoordinated  through  grain  growth,  i.e.,  grain  growth  is  not  required 
for  densification  of  this  special  array. 

In  general,  because  particles  are  not  packed  in  a  symmetrical 
arrangement,  the  coordination  number  for  pores  can  be  widely  distributed. 

This  is  particularly  true  for  real  powders  that  are  packed  as  agglomerates.6 
In  addition,  because  of  their  nonsymmetric  arrangement,  particles  rearrange 
during  the  Initial  stages  of  neck  growth  causing  some  pores  to  increase  their 
coordination  number.4  It  is  suspected  that  a  large  fraction  of  pores  in  a 
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real  powder  compact  have  coordination  numbers  greater  than  critical.  Thus, 
during  initial  sintering,  only  a  fraction  of  the  pores  are  thermodynamically 
unstable  and  disappear  without  grain  growth.  Stable  pores  (i.e.,  those  with 
n  >  nc)  will  shrink  until  grain  growth  can  lower  their  coordination  number  to 
make  them  unstable  and  disappear.  With  grain  growth,  the  stable  pores  will 
shrink  in  proportion  to  the  increased  grain  size  with  the  same  relation  as 
observed  for  the  ring  of  spheres.  If  the  void  volume  is  equally  distributed 
between  pores  with  coordination  numbers  <  8nc,  then  a  grain  size  increase  of 
(2 1/3) 3  lead  tQ  con)p-|ete  densif ication  and  a  near  linear  grain  size  vs 

densif ication  relation  similar  to  that  shown  in  Fig.  9  for  the  ring  array. 

If,  on  the  other  hand,  the  void  volume  was  weighted  toward  smaller  coordina¬ 
tion  numbers,  then  most  of  the  densification  will  only  require  a  modest 
increase  in  grain  size,  while  the  last  and  smallest  increase  in  density  will 
require  the  greatest  increase  in  grain  size.  Thus,  the  extent  of  grain  growth 
required  for  densification  will  depend  on  the  pore  coordination  number  distri¬ 
bution  function.  Assuming  that  pores  remain  coordinated  by  grains  (i.e., 
pores  are  not  swallowed  up  by  abnormally  growing  grains),  the  grain  size 
required  to  achieve  theoretical  density  will  depend  on  the  pore  with  the 
largest  coordination  number. 

With  the  recognition  that  grain  growth  leads  to  densification,  it 
should  also  be  recognized  that  densification  kinetics  will  depend  on  both 
sintering  kinetics  (intraparticle  mass  transport  to  the  neck  region)  and  grain 
growth  kinetics  (interparticle  mass  transport).  If  the  path  length  for  inter¬ 
particle  mass  transport  is  larger  and/or  the  driving  force  for  interparticle 
mass  transport  smaller  than  those  for  sintering,  grain  growth  kinetics  will  be 
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the  rate-limiting  process  and  dominate  the  densification  kinetics.  In  addi¬ 
tion,  grain  growth  will  increase  the  path  lengths  for  both  intra-  and  inter¬ 
particle  mass  transport.  This  increase  in  path  length  is  prescribed  by  the 
amount  of  grain  growth  required  to  reduce  the  coordination  numbers  of  all 
pores  to  their  critical  number  and  is  thus  related  to  the  coordination  number 
distribution  function.  Theoretical  and  experimental  kinetic  studies  of 
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densification  should  be  reevaluated  with  these  thoughts  in  mind. 
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Fig.  1  Configurational  changes  during  sintering  and  interparticle  mass 

transport  leading  to  special  configuration  (d)  where  grain  boundary 
can  move  through  small  grain  and  decrease  the  free  energy  of  the 
system. 
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(a)  Geometrical  parameters  needed  to  define  the  radii  of  the  grains 
and  grain  boundary,  (b)  Geometrical  parameters  that  define  the  spe¬ 
cial  configuration  for  boundary  motion  through  the  smaller  grain. 
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Fig.  4  Configurational  changes  produced  during  sintering  and  interparticle 
mass  transport  for  three  colinear  particles.  Note  that  disappearance 
of  smaller  grain  (configuration  c)  reinitiates  sintering. 
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(a)  (b)  (c) 


(d)  (e) 


Fig.  5  Configurational  changes  produced  during  sintering  and  interparticle 
mass  transport  for  three  noncolinear  particles.  Note  that  sintering 
is  reinitiated  first  at  one  point  (configuration  c). 
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Fig.  6  Geometrical  parameters  needed  to  define  contact  angle  where  two 
larger  grains  touch  during  sintering  of  noncolinear  particles. 
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Shrinkage  strain  for  linear  array  of  two  different  sets  of  particles 
as  a  function  of  dihedral  angle.  Lower  curve  represents  shrinkage 
when  particles  sinter.  Middle  curve  represents  shrinkage  when  every 
other  particle  disappears  by  interparticle  mass  transport  and  bound¬ 
ary  motion.  Upper  curve  is  shrinkage-produced  when  remaining  grains 


SC5410.FR 


THERMODYNAMICS  OF  DE NS  I F I  CAT  ION ,  PART  III: 
EXPERIMENTAL  RELATION  BETWEEN  GRAIN  GROWTH  AND  PORE  CLOSURE 


Bruce  Kellett 

Department  of  Materials  Science  and  Engineering 
University  of  California  -  Los  Angeles 
Los  Angeles,  California  90024 

F.F.  Lange 
Structural  Ceramics 
Rockwell  International  Science  Center 
1049  Camino  Dos  Rios 
Thousand  Oaks,  California  91360 


Rockwell  International 

Science  Center 


THERMODYNAMICS  OF  DENS  IF ICAT ION ,  PART  III: 
EXPERIMENTAL  RELATION  BETWEEN  GRAIN  GROWTH  AND  PORE  CLOSURE 


Bruce  Kellett 

Department  of  Materials  Science  and  Engineering 
University  of  California  -  Los  Angeles 
Los  Angeles,  California  90024 

F.F.  Lange 
Structural  Ceramics 
Rockwell  International  Science  Center 
1049  Camino  Dos  Rios 
Thousand  Oaks,  California  91360 


ABSTRACT 

The  effect  of  grain  growth  on  conditions  for  pore  disappearance  was 
studied,  using  an  A^Oj/ZrC^  composite  containing  monosize  pores  introduced  by 
incorporating  plastic  spheres  (1  urn,  2  um  and  4  ym)  into  the  powder  prior  to 
sintering  at  1600°C.  Pores  uncovered  within  several  microns  of  the  surface  by 
careful  polishing  were  characterized  with  respect  to  their  size  and  number 
density  after  different  periods  of  heat  treatment  at  1600°C.  Grain  growth, 
which  occurred  during  heat  treatment,  was  also  characterized.  Results  show 
that  the  pores  maintain  their  size  while  decreasing  their  number  density  dur¬ 
ing  grain  growth.  Pores  that  remained  after  heat  treatment  were  coordinated 
by  more  than  six  grains  in  the  plane  of  polish.  These  results  are  consistent 
with  theoretical  work,  which  suggests  that  pores  are  stable  when  coordinated 
by  more  than  a  critical  number  of  grains.  The  current  experimental  work  sug¬ 
gests  that  pores  quickly  disappear  once  grain  growth  reduces  their  coordina¬ 
tion  number  to  a  critical  value. 
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INTRODUCTION 


Most  sintering  theories  have  been  based  on  ideas  formulated  by 
Herring, 1  who  pointed  out  that  mass  will  be  transferred  to  surfaces  of  lower 
curvature.  Kuczynki's  analysis  of  spherical  pores,  based  on  these  ideas, 
concluded  that  pores  continuously  decrease  in  size  at  a  rate  dependent  on  the 
mode  of  mass  transport.  In  a  polycrystal,  pores  are  only  spherical  if  they 
reside  within  a  single  grain  that  exhibits  isotropic  surface  energy.  Pores 
intersected  by  grain  boundaries  exhibit  surface  curvatures  that  depend  on  the 
number  of  intersecting  grain  boundaries  and  the  dihedral  angle  ( i^e ,  the  dihed¬ 
ral  angle  depends  on  the  ratio  of  surface  and  grain  boundary  energies).  Pre¬ 
vious  sintering  theories  based  on  mass  transport  to  the  contact  region  (region 
of  negative  surface  curvature)  between  two  spherical  particles  do  not  examine 
the  conditions  for  pore  closure  because  two  touching  particles  cannot  define  a 
pore. 

■5 

Kingery  and  Francois  were  the  first  to  consider  the  geometry  of  pore 
closure  in  polycrystalline  materials.  They  suggested  that  the  flow  of  mass 
either  into  or  away  from  a  pore  will  depend  on  its  surface  curvature  and  that 
the  surface  curvature  of  the  pore  would  depend  on  both  the  number  of  grains 
coordinating  the  pore  and  the  dihedral  angle.  Using  geometrical  arguments, 
they  showed  that  if  the  dihedral  angle  was  to  be  maintained  at  each  grain¬ 
boundary/pore-surface  intersection,  the  curvature  of  the  pore  surface  between 
these  intersections  would  depend  on  the  number  of  intersecting  grain  bound¬ 
aries.  They  concluded  that  if  the  number  of  intersecting  boundaries  (number 
of  grains  coordinating  the  pore)  is  less  than  a  critical  value,  then  the 
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pore’s  surface  will  be  concave  (looking  from  within  the  pore),  which  would 
promote  mass  transport  to  the  pore's  surface  and  pore  closure.  They  also  con¬ 
cluded  that  if  the  pore's  coordination  number  is  greater  than  a  critical  val¬ 
ue,  its  surface  will  be  convex,  promoting  mass  transport  away  from  the  pore 
and  thus,  pore  growth.  The  critical  pore  coordination  number  separating  pore 
closure  from  pore  growth  is  dependent  on  the  dihedral  angle.  It  must  be  con¬ 
cluded  that  they  compared  the  pore's  surface  with  an  external,  flat  surface. 

In  a  companion  paper, ^  we  determined  the  lowest  free  energy  config¬ 
uration  for  closed  arrays  of  identical  particles  which  contained  a  single  pore 
(rings  and  polyhedra).  They  showed  that  all  pores  within  closed  arrays  will 
shrink  during  sintering,  but  only  those  with  a  coordination  number  less  than  a 
critical  value  disappear.  Pores  with  a  coordination  number  greater  than  a 
critical  number  will  shrink  to  an  equilibrium  size.  As  shown  in  the  next  sec¬ 
tion,  the  same  conclusion  is  obtained  for  the  case  of  an  isolated  pore  within 
a  bulk  specimen. 

As  shown  in  Fig.  1,  grain  growth  can  significantly  lower  a  pore's 
coordination  number  and  change  its  surface  curvature  from  convex  to  concave. 

It  is  well  known  that  grain  growth  occurs  throughout  all  stages  of  sintering.^ 
It  is  thus  hypothesized  that  grain  growth  can  be  a  desirable  phenomenon  in 
helping  to  eliminate  pores  by  changing  their  mean  surface  curvature.  Experi¬ 
ments  were  designed  to  test  this  hypothesis.  Materials  were  fabricated  con¬ 
taining  spherical  voids,  which  were  the  only  detectable,  remnant  porosity 
after  sintering.  These  materials  were  heat  treated  to  induce  grain  growth. 
Observations  were  made  to  determine  the  conditions  where  the  large,  spherical 
pores  disappeared.  As  will  be  seen,  these  conditions  are  consistent  with  the 
initial  hypothesis. 
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THERMODYNAMIC  ASPECTS  OF  ISOLATED  PORES  RELATIVE  TO  THE  EXTERNAL  SURFACE 

Our  calculations  on  pore  growth  and  shrinkage  are  based  on  a  two- 
dimensional  model  in  which  all  pores,  surrounded  by  identical  grains,  are 
within  a  bulk  specimen  of  radius  (R),  as  shown  in  Fig.  2.  The  dihedral  angle, 
determined  through  Young's  equation  (cos (u^e/2)  =  yb/2ys),  is  identical  at  all 
grain  boundary-surface  junctions. 

As  described  by  Kingery  and  Francois,  pores  with  a  coordination 
number  greater  than  a  critical  value,  i.e.,  n  >  nc,  where 

nc  =  2n/(  n-«pe )  (1) 

have  a  convex  curvature,  pores  with  a  coordination  number  equal  to  the  criti¬ 
cal  value  have  flat  surfaces,  and  pores  with  n  <  nc  have  a  concave  curvature. 
The  external  surface  always  has  concave  curvature. 

Pore  energy  was  equated  to  its  surface  energy  minus  the  grain  bound¬ 
ary  energy  it  removes  from  the  dense  body,  expressed  as 

E  -  YsAs-YbAb  ■  (^) 

The  surface  area  per  unit  thickness  (A$)  and  the  grain  boundary  area  per  unit 
thickness  (Ab)  removed  by  the  pore  are  a  function  of  the  pore  coordination 
(n),  pore  radius  of  curvature  (r  )  and  pore  size  (Rp),  as  shown  in  Fig.  2: 

A$  =  n(2n/n+iiie--n)rp  ,  (3) 

82 

J7003A/sn 


Li  LI1U‘UN.I1.l\J4A^t»i 


Rockwell  International 

Science  Center 


Ab  =  n  Rp  .  (4) 

It  can  also  be  shown  that  pore  volume  per  unit  thickness  (Vp)  is  given  by 

Vp  =  n/2((2ir/n+i|)e-iT+sin(^e+2Ti/n)]rp2  +  si n(2*/n)Rp2)  .  (5) 

Using  the  relationship  between  rp  and  Rp: 

rp  =  Rp  sin(n/n)/cos(ii.g/2+i!/n)  ,  (6) 

the  derivative  of  pore  enegy  with  respect  to  pore  volume  simplifies  to* 

dEpAJVp  =  -Ys/rp  .  (7) 

Equation  (7)  implies  that  when  considered  alone,  isolated  pores  with  convex 
surfaces  decrease  their  energy  by  growing,  and  pores  with  concave  surfaces  de¬ 
crease  their  energy  by  disappearing.**  The  energy  of  pores  with  flat  sides 
(rp  =  »)  is  independent  of  pore  size.  These  conclusions  were  expressed  by 
Kinergy  and  Francois,3  and  taken  to  define  the  condition  for  pore  growth  or 


*Pore  radius  of  curvature  is  positive  for  convex  pore  curvature  and  negative 
for  concave  pore  curvature. 

♦♦Equation  (7)  is  equivalent  to  the  Kelvin  Equation:  u  =  (1/R^  +  l/R?); 

w  chemical  potential,  n  atomic  volume,  R^  and  R2  are  the  orthogonal  radii 
of  surface  curvature. 
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disappearance.  However,  this  cannot  be  correct  because  changes  in  pore  volume 
also  cause  changes  in  the  energy  of  the  external  surface. 

To  calculate  the  change  in  energy  of  the  exterior  surface  with  re¬ 
spect  to  pore  volume,  consider  the  external  region  shown  in  Fig.  2.  The 
change  in  the  energy  of  the  exterior  surface  can  then  be  calculated  as  the  sum 
of  the  change  in  its  surface  and  grain  boundary  energies, 


dEe  = 


dA'  + 
s 


Yb  dAb 


(8) 


where  dA^  is  the  deriva^'ve  of  the  surface  area  per  unit  thickness,  and  dA^  is 
the  derivative  of  the  grain  boundary  area  per  unit  thickness.  It  can  be  shown 
that 


dA^  =  ne(*+2*/ne-^e)dre  ,  (9) 

dA'  =  nedR  ,  (10) 

where  ng  is  the  number  of  grains  surrounding  the  specimen,  and  dr£  is  the 
derivative  of  the  exterior  surface  curvature. 

Using  the  relationship  between  rg  and  R, 

re  =  R  sin(n/ne)/cos(ii»e/2-n/ne)  ,  (11) 

and  Young's  equation,  the  derivative  of  the  exterior  energy  with  respect  to 
the  volume  of  the  specimen  (V$)  simplifies  to: 
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dEe/dV$  =  Ys/re  .  (12) 

Since  the  exterior  surface  curvature  is  positive,  the  exterior  surface  energy 
always  increases  with  increasing  volume. 

If  we  assume  that  the  change  in  pore  volume  is  equal  to  the  change  in 
specimen  volume  (dVp  =  dV$),  then  the  total  change  in  the  system's  energy  with 
respect  to  pore  volume  is  the  sum  of  Eqs.  (7)  and  (12),  and  is  a  function  of 
the  radius  of  curvature  of  the  pore's  surface  (rp),  and  the  radius  of  curva¬ 
ture  of  grains  on  the  external  surface  (r0): 

dEt/dVp  =  Ys(l/re-l/rp)  .  (13) 

For  concave  pores  (rp  <  0),  the  energy  of  the  system  decreases  with  pore 
volume,  implying  that  these  pores  will  continuously  shrink  and  disappear, 
kinetics  permitting.  Convex  pores  (rp  >  0)  will  either  shrink  or  grow  depend¬ 
ing  on  the  relative  values  of  the  curvatures,  i.e.,  when  rg  >  rp  pore  growth 
decreases  system  energy,  and  for  rg  <  rp,  pore  shrinkage  decreases  system 
energy.  Similar  conclusions  were  obtained  by  Cannon6  through  a  different 
approach. 

Equilibrium  (dEt/dVp  =  0)  occurs  when: 

re  =  rp  »  04) 
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suggesting  that  pores  with  convex  surfaces  will  be  stable  when  their  curvature 
equilibrates  with  that  of  the  exterior  surface. 

The  exterior  radius  of  curvature  can  be  equated  to  the  grain  size  (D) 
(see  Fig.  2) : 

Re  =  D/[2cos(ve/2) ]  .  (15) 

Using  Eqs.  (6),  (14)  and  (15),  it  can  be  shown  that  at  equilibrium,  the  pore 
size  is  dependent  on  the  grain  size  (D)  and  the  pore  coordination  number  (n): 


Rp.O 


“s  <r  *  H 

2cos  (jjr)  sin(J) 


(16) 


Equation  (16)  shows  that  when  n  <  nc  =  2ti/ (n-vpe) ,  Rp  =  0  (the  pore  disap¬ 
pears),  and  when  n  >  nc,  the  pore  is  stable  with  a  finite  size  (Rp  >  0). 
Equation  (16)  also  shows  that  if  grain  growth  occurs  without  changing  the 
pore's  coordination  number,  the  pore  will  grow.  On  the  other  hand,  if  grain 
growth  lowers  the  coordination  of  the  stable  pore,  a  condition  will  be  reached 
such  that  the  pore  will  become  unstable  and  will  disappear. 

These  results  differ  from  those  of  Kingery  and  Francois  who  stated 
that  pores  of  n  >  nc  experience  a  driving  force  for  growth. 
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EXPERIMENTAL 

Pores  of  uniform  size  were  introduced  into  a  sintered  body  by  mixing 
a  small  fraction  (<  1  vol%)  of  polystrerene  microspheres  of  narrow  size  dis¬ 
tribution  (~  1%)  into  a  flocced,  two-phase  (A^Oj/ZrC^)  powder  slurry.  The 
colloidal  method  of  treating  the  powders  to  decrease  agglomerate  size  and  the 
method  of  mixing  has  been  reported  elsewhere.^  The  slurry  was  consolidated  by 
slip  casting  and  was  sintered  after  drying  at  1600°C.  The  plastic  micro¬ 
spheres  burned  out  prior  to  the  initiation  of  bulk  shrinkage.  The  resulting 
microstructure  contained  spherical  pores,  diminished  in  size  due  to  bulk 
shrinkage,  highly  coordinated  by  smaller  grains.  Different  specimens  were 
fabricated  containing  pores  produced  with  1  pm,  2  pm  and  4  pm  plastic  spheres. 
Grain  growth  was  induced  by  reheating  to  the  sintering  temperature  for  differ¬ 
ent  periods.  Grain  size  and  pore  size  distributions  and  pore  number  densities 
were  determined  as  a  function  of  heat  treatment  period. 

Preliminary  work  with  single  phase  A1 2O3  powder  produced  abnormal 
grain  growth  at  the  heat  treatment  temperature  causing  many  of  the  pores  to 
relocate  within  the  larger  abnormal  grains,  preventing  observations  on  the 
influence  of  grain  growth  on  pore  closure.  For  this  reason,  10  vol%  Zr02 
(+6.6  mole%  V2O3)  was  incorporated  with  the  AT 203  to  prevent  abnormal  grain 
growth,  as  described  elsewhere.® 

After  each  heat  treatment,  pores  located  several  microns  beneath  the 
surface  were  uncovered,  without  producing  significant  pull-outs,  by  careful 
diamond  polishing  with  sequential  grit  sizes  of  4,  1  and  1/4  pm.  After  the 
specimens  were  polished,  they  were  thermally  etched  at  1450°C  to  reveal  grain 
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boundaries.  Scanning  electron  micrographs  were  taken  at  a  magnification  suit¬ 
able  to  obtain  statistical  information  concerning  grain  size,  pore  size  and  to 
determine  the  number  of  pores  per  unit  area  (pore  number  density).  The  appro¬ 
priate  microstructural  features  were  traced  separately  onto  a  transparent 
overlay  to  allow  the  use  of  an  image  analyzer*  without  the  ambiguous  defini¬ 
tion  of  grain  and  pore  boundaries.  The  software  associated  with  the  image 
analyzer  converted  the  areas  of  the  microstructural  features  to  an  equivalent 
circle;  grain  size  and  pore  size  are  reported  as  equivalent  circle  diameters. 

RESULTS 


Figure  3  illustrates  that  the  median  grain  size  of  both  the  A1203  and 
Zr02  grains  in  the  Al203/Zr02  composite  was  observed  to  increase  linearly  with 
respect  to  log  time  for  heat  treatment  at  1600°C.  The  smaller  Zr02  inclu¬ 
sions,  primarily  located  at  four-grain  junctions,  grew  by  coalescence;^  the 
ratio  of  the  median  sizes  of  the  two  phases  was  constant  (~  0.6)  throughout 
the  heat  treatment.  The  three  grain  size  histograms  shown  in  Fig.  4  illus¬ 
trate  that  the  increased  median  grain  size  is  due  principally  to  the  broaden¬ 
ing  of  the  grain  size  distribution  with  increasing  heat  treatment  periods. 

Figure  5  illustrates  two  pore  size  histograms  (heat  treatment  periods 
of  1  and  13  h)  for  material  fabricated  with  the  2  urn  plastic  spheres.  The 
smaller  median  size  of  the  pores  relative  to  the  size  of  the  starting  plastic 
spheres  reflects  the  linear  shrinkage  of  the  powder  compact  (i.e.  all  pores 
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shrink  during  sintering)  and  the  statistical  mean  diameter  intersected  by  the 
polished  plane.  As  shown,  the  pore  size  distribution  does  not  significantly 
change  during  heat  treatment. 

The  median  pore  size  and  the  pore  number  density  plotted  as  a  func¬ 
tion  of  grain  size  are  illustrated  in  fig.  6a  and  6b  for  the  materials  fabri¬ 
cated  with  the  three  different  sized  plastic  spheres.  As  shown,  the  median 
pore  size  remains  constant,  wnereas  the  pore  number  density  decreases  with 
increasing  grain  size.  The  decrease  in  pore  number  density  for  the  larger 
pores  produced  with  4  um  plastic  spheres  was  less  dramatic  than  the  smaller 
pores.  The  slight  increase  in  pore  size  at  the  largest  grain  size  observed 
for  the  smaller  pores  (produced  with  the  1  pm  and  2  urn  spheres)  is  not  con¬ 
sidered  significant  since  as  the  pore  coordination  number  decreased,  fewer 
pores  were  observed  and  their  shape  became  more  irregular,  making  them  more 
difficult  to  distinguish  from  the  occasional  pull-out  of  Zr02  inclusions 
during  polishing. 


DISCUSSION 

Pores,  remnants  of  plastic  microspheres,  were  observed  to  remain  at  a 
uniform  size  while  their  number  density  decreased  with  increasing  grain  size. 
This  observation  implies  that  pores  do  not  continuously  shrink,  but  are 
quickly  removed  in  a  discontinuous  manner.  The  thermodynamics  of  pore  stabil¬ 
ity  suggests  that  pores  coordinated  by  a  greater  than  critical  number  of 
grains  are  stable  and  can  only  disappear  when  their  coordination  number  be¬ 
comes  equal  to  the  critical  number.  Since  the  mass  transport  distance  between 
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the  observed  pore  surface  and  the  external  surface  was  no  greater  than  several 
pore  diameters  in  the  experiments  described  above,  once  the  critical  coordi¬ 
nation  condition  was  achieved  by  grain  growth,  pore  disappearance  could  be 
rapid.  Because  real  polycrystall ine  materials  exhibit  a  grain  size  distribu¬ 
tion,  the  coordination  of  some  pores  will  decrease  to  the  critical  number  be¬ 
fore  others;  thus,  pore  number  density  would  be  expected  to  decrease  grad¬ 
ually.  Remaining  pores  would  have  a  coordination  number  greater  than  the 
critical  value. 

Figure  7  illustrates  the  coordination  number  of  the  pores  remaining 
after  heat  treatment,  determined  on  the  polished  plane,  as  a  function  of  grain 
size  for  materials  fabricated  from  three  different  sphere  sizes.  The  coordi¬ 
nation  number  of  the  pores  produced  with  the  2  pm  and  4  pm  spheres  initially 
decreases  with  increasing  grain  size,  whereas  the  coordination  number  of  the 
smaller  pores  produced  with  the  1  pm  spheres  remains  nearly  constant  through¬ 
out  grain  growth.  For  the  case  of  the  pores  produced  from  the  1  pm  and  2  pm 
spheres,  the  smallest,  average  coordination  number  achieved  at  larger  grain 
sizes  was  6,  whereas  that  for  the  larger  pores  produced  with  the  4  pm  spheres 
was  slightly  larger. 


-  j 
■*  *  * 


) 
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In  an  approximate  way,  one  could  relate  this  in-plane  coordination  to 
the  coordination  of  a  dissected  polyhedron;  i.e.,  an  in-plane  coordination  of 
5  or  6  would  result  from  the  bisection  of  a  dodecahedron  (pore  coordination 
12),  and  from  this  polyhedron  predict  the  equilibrium  dihedral  angle  (for  a 
dodecahedron  shaped  pore  the  dihedral  angle  is  138°).  However,  grains  are  not 
synuietrically  arranged  around  real  pores.  Also,  the  dihedral  angle  would  be 
expected,  for  an  anisotropic  material,  to  depend  on  the  relative  orientations 
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of  adjoining  grains.  Thus,  it  would  be  difficult  to  determine  the  critical 
coordination  number  with  the  current  data. 

A  reasonable  conclusion  is  that  if  pores  within  the  initial  powder 
compact  have  a  coordination  number  greater  than  the  critical  value,  then  grain 
growth  would  be  a  desirable  phenomenon,  since  it  causes  stable  pores  to  dis¬ 
appear  by  decreasing  their  coordination  number.  Thus,  the  initial  distribu¬ 
tion  of  pore  coordination  numbers  within  a  crystalline  powder  compact  will 
determine  the  amount  of  grain  growth  required  to  achieve  a  dense  material. 
Pores  that  can  not  be  eliminated  through  grain  growth  must  be  removed  by  de¬ 
formation,  e.g.  hot-pressing  or  HIPing. 
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Surface  curvature  for  two  pores  of  coordination  (n)  8  and  3,  3nd  a 
region  of  the  exterior  surface. 
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Pore  size  histogram  of  the  material  fabricated  from  2  urn  plastic 
spheres  after  heat  treatments  of  a)  1  and  b)  13  h. 
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1  um  MODEL  PORES 


2  yum  MODEL  PORES 


4  nm  MODEL  PORES 


MEDIAN  GRAIN  SIZE 

(a)  Median  pore  size  and  (b)  number  pore  density  of  material  fabri¬ 
cated  from  1,  2,  and  4  um  plastic  spheres  as  a  function  of  the  com¬ 
bined  AloOo/ZrOp  median  grain  size,  for  heat  treatments  at  1600°C. 
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MEDIAN  GRAIN  SIZE 


Pore  coordination,  in  the  plane  of  polish,  of  the  1,  2,  and  4  um 
plastic  spheres  fabricated  material,  as  a  function  of  the  combined 
A^Oj/ZrC^  medium  grain  size,  for  heat  treatments  at  1600°C. 
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ABSTRACT 

The  effect  of  stress  on  pore  closure  was  studied  by  HIPing  (Hot 
isostatic  Pressing)  and  hot  forging  sintered  A 1 2O3 ,  cubic-Z^,  and  finer 
grained  A^Oj/ZK^  (cubic)  composite  materials  containing  pores  remnant  of 
monosized  plastic  spheres.  For  the  temperature-time  conditions  explored, 
plastic  deformation  was  the  dominate  mechanism  for  pore  shape  change  and 
closure.  Namely,  only  pores  in  the  highly  deformable  two-phase  material  where 
observed  to  either  change  shape  (under  uniaxial  forging)  or  uniformly  shrink 
(during  HIPing).  Large  plastic  strains  (0.3  to  0.6,  dependent  on  location) 
were  required  to  produce  pore  closure  under  forging  conditions.  At  tempera¬ 
tures  where  the  two-phase  materials  were  sufficiently  deformable  to  affect 
pore  closure,  pores  near  the  surface  would  puncture  as  the  thin  web  of 
material  separating  the  pore  from  the  surface  was  deformed  into  the  pore. 
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These  results  are  discussed  in  terms  of  the  common  practice  of  HIP  finishing 
sintered  materials. 


1.0  INTRODUCTION 


n 


Even  though  HIPing  and  hot  forging  has  long  been  recognized  as  a 
technique  for  increasing  densif ication  kinetics,  and  end  point  density,  disa¬ 
greement  exists  on  the  mechanism(s)  causing  pore  shrinkage  and  disappearance. 
Two  different  viewpoints  have  been  used  to  describe  the  removal  of  isolated 
pores,  namely  stress  assisted  sintering  and  creep  induced  pore  collapse. 

Coble1  assumed  that  isolated  pores  within  a  polycrystalline  material 
could  be  modeled  as  spherical  voids  of  radius  rp,  and  reasoned  that  an  applied 
compressive  stress  (oa)  would  decrease  the  concentration  of  vacancies  at  grain 
boundaries  and  hence  increase  diffusional  transport  of  mass  to  the  void  to  aid 
in  pore  closure.  He  concluded  that  the  applied  stress  would  enhance  the 
driving  "stress"  (oj)  for  pore  closure  as  expressed  by 


°d  =  °a  '  Vrp 


where  ys  is  the  surface  energy  of  the  pore  per  unit  area.  Since  Coble  assumed 
that  pores  always  have  a  negative  surface  curvature  (i.e.,  spherical  voids) 
all  pores  were  thermodynamically  unstable  and  would  disappear.  An  applied 
stress  would  simply  aid  in  their  disappearance. 
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But,  as  pointed  out  by  Kingery  and  Francois/  the  surface  curvature 
of  pores  within  a  polycrystali ine  material  depends  on  the  number  of  coordinat¬ 
ing  grains  and  the  ratio  of  the  grain  boundary  to  surface  energies  as 
expressed  through  Young's  equation.  They  showed  that  when  the  number  of 
coordinating  grains  is  greater  than  a  critical  number  (n  >  nc)  the  pore's 
surface  is  convex  (rp  >  0),  and  when  n  <  nc,  the  surface  is  concave  (looking 
from  within  the  pore)  and  rp  <  0.  When  the  free  energy  of  a  body  containing 
an  isolated  pore  is  determined  as  a  function  of  the  pore  volume,  it  can  be 

O 

shownJ  that  when  n  >  nc,  the  isolated  pore  will  either  shrink  or  grow  until 
its  curvature  is  identical  to  the  surface  curvature  of  grains  on  the  external 
surface.  That  is,  only  pores  with  n  <  nc  (rp  <  0)  are  thermodynamically 
unstable  and  disappear,  kinetics  permitting.  With  this  in  mind,  it  can  be 
shown  (see  Appendix)  that  the  effect  of  an  applied  stress  on  the  driving 
"stress"  for  pore  closure  can  be  expressed  as 


°d  =  2ys 


(i-  -  ±-) 
\r  r  ) 


+  0 . 


(2) 


where  re  is  the  radius  of  curvature  of  the  grains  on  the  external  surface. 
Since  rp  changes  with  pore  volume  (the  absolute  magnitude  of  rp  decrease  with 
pore  volume),  whereas  re  only  depends  on  grain  size,  it  can  be  shown  (see 
Appendix)  that  the  equilibrium  size  Rp  of  the  isolated  pore  changes  with 
applied  pressure  as 
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R_  = 


cos  +  J) 


°a  .  2  *e. 


(3) 


^  +  D  cos  r)  Sin  n 


Equation  (3)  shows  that  an  applied  pressure  will  simply  decrease  the  equilib¬ 
rium  size  of  the  stable  pore  (n  >  nc),  whereas  unstable  pores  (n  <  nc)  are 
expected  to  disappear  faster  with  an  applied  pressure. 

Deformation  induced  pore  collapse  is  chiefly  dependent  on  the 
material’s  ability  to  deform  in  a  plastic  manner  and,  in  particular,  on  the 
material's  stress-strain  rate  behavior.  Wilkinson  and  Ashby,4  assuming  power 
law  creep  behavior  (i.e.  o  =  Atn),  determined  that  spherical  pores  under  pure 
hydrostatic  conditions  remain  spherical  and  collapse  at  a  constant  volume 
rate.  More  recent  work  by  Budiansky,  et  a!,5  determined  that  under  forging 
conditions  pores  collapse  into  penny  shaped  cracks,  as  one  might  expect.  It 
was  also  found  that  slight  perturbations  from  hydrostatic  stress  also  cause 
initially  spherical  pores  to  collapse  into  either  a  penny  (crack)  or  rod 
shaped  pores.  Perturbations  from  hydrostatic  stress  might  arise  during  HIPing 
(e.g.,  from  the  interaction  of  the  stress  state  about  neighboring  pores), 
suggesting  that  pore  closure  during  HIPing  might  proceed  by  a  pore  shape 
change  instability. 

To  determine  the  effects  of  stress  and  strain  on  pore  closure, 
sintered  specimens,  fabricated  with  mono-sized  pores  with  n  >  nc,  were  both 
HIPed  at  different  temperatures,  and  hot  forged  to  different  strains.  Micro¬ 
scopic  examination  as  used  to  follow  pore  closure  conditions.  Results 
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indicate  that  conditions  that  promote  bulk  deformation  are  much  more  effective 
for  pore  closure  relative  to  stress  enhanced  sintering. 


2.0  EXPERIMENTAL  PROCEDURES 


Pores  of  uniform  size  were  created  by  mixing  a  small  fraction  of 
plastic  spheres,  manufactured  to  a  narrow  size  distribution,*  into  a  ceramic 
powder  slurry  which  was  colloidally  processed  to  remove  large  agglomerates  as 
described  elsewhere.®  Both  single-phase  AI2O3**  and  ZrOpt  and  a  two-phase 
mixture  of  these  Al203/Zr02  powders  were  used  (the  Zr02  powder  was  prealloyed 
with  12  weight%  (6.6  mole%)  Y2O3,  which  produces  the  cubic  structure  after 
densif ication) .  The  slurry  was  ultrasonicated  in  the  dispersed  state,  to 
insure  good  mixing,  and  then  flocced  to  prevent  separation  (due  to  differences 
in  the  settling  rates  between  the  plastic  spheres  and  the  ceramic  powder). 
Slurries  were  consolidated  by  pressure  filtration,  as  described  elsewhere,7 
dryied  overnight  in  an  oven,  and  sintered  in  air  for  one  hour  at  either 
1600°C,  for  the  A^C^/ZrC^  composites  and  A1 203  material,  and  1500°C  for  the 
Zr02  material.  The  plastic  spheres  burn  out  prior  to  bulk  shrinkage,  leaving 
highly  coordinated  pores.  Different  specimens  containing  different  size 
pores,  remnant  of  the  different  size  plastic  spheres  (4  to  50  ym)  where 
fabricated.  Specimens  were  diamond  cut  and  polished. 


*Duke  Scientific,  445  Sherman  Ave.,  Palo  Alto,  CA  94306, 
♦♦Alcoa  A  16-SG. 

JZircar  Products  Inc.,  110  North  Main  St.,  Florida,  NY  10921. 
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HIPing  was  performed  in  a  MiniHIP*.  The  pressure-temperature 
schedule  consisted  of  an  initial  pressurization  to  70  MPa,  ramping  to 
temperature  (either  1400°C,  1450°C,  1475°C,  and  1500°C),  a  temperature  hold 
for  10  min,  followed  by  cooling  and  depressurization  to  ambient  conditions. 

Hot  forging  was  performed  with  a  cantilevered  dead-weight  creep 

O 

apparatus,  previously  described.  Uniaxial  stresses  were  applied  to  the 
specimen  with  SiC  rams  through  polished  silicon  carbide  pads.  Forging  was 
done  at  constant  temperature  (1500‘C)  and  at  constant  load.  The  load  bearing 
area  of  the  specimen  increases  during  forging  decreasing  the  actual  stress  on 
the  specimen.  Also,  frictional  constraint  between  the  specimen  and  the 
loading  pads  create  hydrostatic  stresses  within  the  specimen  which  diminishes 
toward  the  specimen  free  surface  (friction  hill  stress).9  The  lack  of 
significant  specimen  barreling  suggests  that  friction  hill  stresses  were 
relatively  minor. 

Forged  or  hiped  specimens  were  diamond  cut  and  polished  for  micro- 
structural  observations. 


*Conaway  (now  manufactured  by  ASEA,  Sweden). 

106 

J7964TC/bw 


1  *  m ~  m 


.vV) 

1 


V/ 


TOW  A  3  H  W  AWMUUUU  U  *? 1 

s 

t 


F. ■,«  PI f.i  ^7  ■.  ■  ■  '.«  ■-■■.■  .'  V  7','  ’? T7  1 


"V  ~_>  ' 


fl* 


Rockwell  International 

Science  Center 


’>  •:  y/ 

v.-s/ 


7.V..' 


( 


3.0  RESULTS 

3. 1  Deformation  Characteristics 

Figure  1  illustrates  the  deformation  behavior  of  the  two  single  phase 
materials  and  one  of  the  two-phase  materials  for  data  obtained  under  hot- 
forging  conditions  at  1500°C  under  a  constant  load  that  produced  an  initial 
stress  of  70  MPa.  Both  single  phase  materials,  which  had  a  larger  average 
grain  size  (3  pm,  A1203;  7.5  pm,  Zr02)  failed  in  ternary  creep,  whereas  the 
fine  grained  (1  pm)  two-phase  materials  exhibited  extensive  deformation 
without  developing  tears  or  cracks.  The  deformation  behavior  of  the  two-phase 
materials  is  detailed  elsewhere.10 

3.2  HIPinq  Experiments 

Optical  micrographs  of  diamond  cut  and  polished  surfaces  of  HIPed, 
Al203/20  vol%  Zr02  specimens  containing  pores  remnant  of  the  25  pm  plastic 
spheres  are  shown  in  Fig.  2.  These  micrographs  show  that  the  pore  number 
density  and  pore  size  decrease  with  increasing  temperature.  The  pores  within 
the  single  phase  materials  showed  no  measurable  change  over  this  same 
temperature  range. 

Figure  3  shows  a  polished  surface  of  a  Al203/20  voU  Zr02  specimen 
(a)  before  and  (b)  after  HIPing  and  illustrates  that  sub-surface  pores 
collapse  to  form  craters.  Normarsky  phase  interference  microscopy  illustrates 
that  the  craters  are  much  larger  than  the  pores.  As  expected  pores  previously 
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open  to  the  surface  do  not  close.  SEM  micrographs  of  these  craters  (Fig.  4) 
Illustrate  that  while  some  craters  are  a  local  depression  in  the  surface, 
others  have  punctured. 

The  single  phase  materials  did  not  exhibit  cratering. 

Occasional  pores  appreared  slightly  elliptical  in  shape  suggestive  of 

c 

the  perturbation  suggested  by  Budiansky  et  al. 

3.3  Hot  Forging  Experiments 

SEM  micrographs  of  polished  and  thermally  etched  surfaces  (Fig.  5) 
show  pores  remnant  of  4  urn  plastic  spheres  located  near  the  edge  of  a  two 
phase  specimen  (containing  10  vol%  Zr02)  where  triaxial  stresses  generated  by 
frictional  constraint  at  the  loading  pads  are  expected  to  be  insignificant. 
Pores  develop  a  crack  like  morphology  with  increased  applied  strain,  and 
significantly,  these  relatively  small  pores  are  still  present  after  an  applied 
engineering  strain  of  0.3.  These  large  strains  were  not  obtainable  for  the 
single  phase  materials. 

Larger  pores,  remnant  of  50  urn  plastic  spheres,  show  (Fig.  6)  similar 
behavior  in  more  detail.  These  pores,  located  near  the  edge  of  the  specimen, 
maintain  a  uniform  width  normal  to  the  applied  load  while  their  height  in  the 
direction  of  the  applied  load  decreases.  Contact  between  the  top  and  bottom 
of  the  pore  is  followed  by  complete  pore  disappearance  for  an  applied 
engineering  strain  of  approximately  0.6. 
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Pores  located  toward  the  center  of  the  forged  specimen,  show  similar 
morphological  changes,  but  appear  to  disappear  at  somewhat  lower  strains 
(about  0.3)  It  is  expected  that  pores  close  to  the  center  of  the  specimen 
might  experience  a  hydrostatic  stress  (generated  by  frictional  constraint)  in 
addition  to  the  applied  uniaxial  stress. 

Pores  in  a  two-phase,  Zr02/20  vol%  2O3 •  material,  which  also 

exhibits  high  ductility  at  1500°C  ,  albeit  somewhat  less  than  the  alumina  rich 
compositions,  developed  cracks  where  tensile  stresses  are  expected  to  develop 
during  forging.  Smaller  tears  and  cavities,  developed  along  the  top  and 
bottom  of  the  pores  (location  of  tensile  stresses)  are  commonly  observed  in 
the  more  deformable  materials  (see  Fig.  6). 

4.0  CONCLUSIONS 


Pore  disappearance  only  occurred  in  the  fine  grain  alumina-zirconia 
material  which  exhibited  extensive  deformation  at  the  hiping  and  forging  temp¬ 
erature.  Cratering  and  puncturing  of  near  surface  pores  during  HIPing  and  the 
observed  change  in  pore  shape  during  forging  strongly  suggests  that  plastic 
deformation  was  responsible  for  pore  closure. 

Pore  closure  during  forging  appears  to  occur  at  an  engineering  strain 
of  approximately  0.6  near  the  free  surface,  and  0.3  near  the  center  of  the 
specimen  where  hydrostatic  stresses  generated  by  frictional  constraint  are 
expected  to  be  greater.  These  results  are  consistent  with  those  of  Budiansky 
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et  al^  who  determined  that  pore  closure  should  occur  at  an  engineering  strain 
of  0.46  under  uniaxial  conditions,  decreasing  to  0  for  hydrostatic  condition. 

The  current  observations  strongly  suggest  that  HIP  finishing  of 

sintered  materials,  commonly  used  to  decrease  the  flaw  populations  associated 

11  10 

with  irregular  voids,  ’  should  be  carried  out  in  a  pressure-temperature 
regime  where  the  material  exhibits  extensive  plastic  deformation.  In  addi¬ 
tion,  although  it  is  commonly  recognized  that  pores  at  the  surface  will  not  be 
closed  during  HIPing,  the  current  observations  show  that  near  surface  pores 
can  create  extensive  surface  damage  during  HIPing. 
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APPENDIX 

EFFECT  OF  AN  EXTERNAL  PRESSURE  ON  THE  EQUILIBRIUM  BETWEEN 
A  PORE  AND  AN  EXTERNAL  SURFACE 

■3 

Previous  calculations  concerning  the  free  energy  of  a  polycrystal¬ 
line  body  of  uniform  grain  size,  D,  containing  an  isolated  pore  showed  that 
the  derivative  of  the  total  free  energy  of  the  system  (pore  plus  external 
surface)  with  respect  to  the  volume  of  the  pore  (Vp)  can  be  expressed  as 

(JE/dVp  =  vs(ke  -  kp)  ,  (Al) 

where  ke  and  kp  are  the  surface  curvatures  of  the  external  grains  and  the 
grains  surrounding  the  pore,  respectively  and  ys  is  the  energy  per  unit  area 
of  the  surface.  ke  is  always  positive.  It  can  also  be  shown  that  the  abso¬ 
lute  magnitude  of  kp  decreases  with  the  size  of  the  pore  (Rp).  It  was  shown 
that  when  the  number  of  grains  coordinating  the  pore  is  greater  than  a  criti¬ 
cal  number  (e.g.,  n  >  nc),  kp  is  also  positive.  Thus,  when  the  coordination 
number  of  an  isolated  pore  is  greater  than  the  critical  value,  the  pore  will 
either  shrink  or  grow  until  the  surface  curvature  of  the  surrounding  grains  is 
identical  to  that  of  the  surface  grains  (viz.  kp  =  ke).  The  surface  curvature 
of  grains  surrounding  pores  with  a  coordination  number  less  or  equal  to  the 
critical  number  (n  <  nc)  is  negative.  Such  pores  are  always  unstable  and 
disappear,  kinetics  permitting.  For  symmetric  pores,  nc  =  2h/(h  -  ii>e),  where 
i|>e  is  the  dihedral  angle  which  defines  the  ratio  of  the  grain  boundary  and 
surface  energies  through  Young's  equation. 
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The  effect  of  an  applied  pressure,  Pa,  on  pore  stability  can  be 
ascertained  by  noting  that  dE/dVp  =  Pa.  Combining  this  term  with  Eq.  (Al) 
results  in 


dE/dVp  =  Y$(ke  -  kp)  ♦  Pa  .  (A2) 

Using  the  relation  between  the  radius  of  curvature  for  the  grains  on  the 
external  surface  and  grain  size  (e.g.,  re  =  D 1 2  cos( «Pe/2) J and  the  relation 
between  pore  radius  (Rp)  and  pore  volume  (see  Ref.  3),  one  can  obtain  the 
following  expression  between  equilibrium  pore  radius  ('i.e.,  when 
dE/dVp  =  0)  and  the  applied  pressure 


(A3) 


Equation  (A3)  shows  that  when  (n  >  nc),  the  pore  always  has  a  finite  size. 


112 

J7964TC/bw 


* 


Rockwell  International 

Science  Center 


1. 


REFERENCES 


R.L.  Coble,  "Diffusion  Models  for  Hot  Pressing  with  Surface  Energy  and 
Pressure  Effects  as  Driving  Forces,"  J.  Appl.  Phys.  £1  | 12 1  4798  (1970), 


*  V  v  « 

■'Nv*. 

V/  v* 

TO, 

i .  .1  ■*> 


2.  W.D.  Kingery  and  B.  Francois,  "Singering  of  Crystalline  Oxides,  I. 
Interaction  Between  Grain  Boundaries  and  Pores,"  pp  471-98  in  Sintering 
and  Related  Phenomena.  Ed.  by  G.C.  Kuczynke,  N.A.  Hooton,  and  G.F. 

Gibbon,  Gordon  Breach,  New  York  (1967). 

3.  B.J.  Kellett  and  F.F.  Lange,  "Thermodynamics  of  Densification,  Part  III: 
Experimental  Relation  Between  Grain  Growth  and  Pore  Closure,"  submitted 
to  the  J.  Am.  Ceram.  Soc. 

4.  D.S.  Wilkinson  and  M.F.  Ashby,  "Pressure  Sintering  by  Power  Law  Creep," 
Act^  Met.  23,  1277  (1975). 

5.  B.  Budiansky,  J.W.  Hutchinson,  and  S.  Slutsky,  "Void  Growth  and  Col  lap-,, 
in  Viscous  Solids,"  in  Mechanics  of  Solids:  The  Rodney  Hill  60th 
Anniversary  Volume,  Ed.  by  H.G.  Hopkins  and  M.J.  Sewell,  Pergamon  Press, 
Oxford. 

6.  F.F.  Lange,  B.I.  Davis,  and  E.  Wright,  "Processing-Related  Fracture 
Origins:  IV.  Elimination  of  Voids  Produced  by  Organic  Inclusions,"  J. 
Am.  Ceram.  Soc.  69  1 1]  66-9  (1986). 

7.  F.F.  Lange  and  K.T.  Miller,  "Kinetics  and  Mechanics  of  Pressure 
Filtration,"  submitted  to  the  J.  Am.  Ceram.  Soc. 

8.  F.F.  Lange,  B.I.  Davis  and  D.R.  Clarke,  "Compressive  Creep  of  SioN4/MgO 
Alloys:  Part  I,"  J.  Mat.  Sci.  15,  600-10  (1979). 

9.  K.R.  Vendatachari  and  R.  Raj,  "Superplastic  Flow  in  Fine-Grain  Alumina," 
J.  Am.  Ceram.  Soc.  69  (2],  135  (1985). 

10.  B.  J.  Kellett  and  F.  F.  Lange,  "Hot  Forging  Characteristics  of 
Transformation  Toughened  A^Oj/Z^  Composites,"  to  be  published. 

11.  V.  Engle  and  H.  Hubner,  "Strength  Improvement  of  Cemented  Carbides  by  Hot 
Isostatic  Pressing,"  J.  Mater.  Sci.,  1_3  [9]  2003-13  (1978). 

12.  F.F.  Lange,  "Processing  Related  Fracture  Origins:  I.  Observations  in 
Sintered  and  Isostactically  Hot-Pressed  Composites,"  J.  Am.  Ceram.  Soc. 

66  [61,  398-8  (1983). 


113 

J7964TC/bw 


»  «  «  -  •«  . 

B 

«  •  • 


•j  -  . 

r:.; 

■>>> 


vw 


s .  .  „ .  >  • 


(nation  strain  (engineering 
obtained  during  forging  ex 
produced  an  initial  stress 


SC 36799 


1 400  °C 


1450°C 


100  fim 


1475°C 


1500°C 


Optical  micrographs  of  cut  and  polished  surfaces  of  A1 2^3/ ( 20  vol  ) 
Zr02  materials  containing  pores  remnant  of  25  pm  plastic  spheres 
sintered  at  1600°C  and  HIPed  at  either  1400°C,  1450°C,  1475°C,  or 
1500°C  at  70  MPa  for  10  min. 
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Fig.  4  Scanning  electron  micrographs  of  HIPed  surface  showing  surface 
dimpling  and  punctures  associated  with  near  surface  pores. 
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Fig.  5  Subterraneous  pores  remnant  of  4  m  plastic  spheres  exposed  by 
polishing  after  different  creep  strains.  Material:  A1 2O3/ 10  vol"' 
Zr02  sintered  at  1600°C  and  forged  at  1500°C. 
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